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Abstract 

^ . Progress on the potential method, recently proposed to investigate hadron interactions in lattice 

QCD, is reviewed. The strategy to extract the potential in lattice QCD is explained in detail. The 
method is applied to extract A^A^ potentials, hyperon potentials and the meson-baryon potentials. 
V^ . A theoretical investigation is made to understand the origin of the repulsive core using the operator 

'sj" I product expansion. Some recent extensions of the method are also discussed. 
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1 Introduction 

Quarks and gluons, the fundamental degrees of freedom of Quantum Chromodynamics (QCD) never 
appear in Nature as asymptotic one-particle states. This phenomena is called conflnement. Instead, 
only their bound states, i.e., mesons and baryons, are observed in experiment and these hadronic states 
are the asymptotic states whose interactions can be parametrized by the S-matrix. An ultimate goal 
for theoretical studies of the strong interaction, therefore, is to extract the properties of the hadronic 
S-matrix from QCD. 

For a description of hadronic interactions, the nuclear force is one of the most fundamental quantities 
in nuclear physics. The origin of the nuclear force, however, is still one of the major unsolved problems 
in strong interaction physics even after the establishment of QCD. Although the nuclear force is not well 
understood theoretically, a large number of proton-proton and neutron-proton scattering data as well 
as deuteron properties have been experimentally accumulated and summarized as the nucleon-nucleon 
(NN) potential, which is designed to reproduce these experimental properties through the Schrodinger 
equation for the NN wave function. Once the potential is determined, it can be used to study systems 
with more than 2 nucleons by using various many-body techniques. 

Phenomenological A^A^ potentials which can flt the A^A^ data precisely (e.g. more than 2000 data 
with x^/dof ~ 1 ) at Tiab < 300 MeV are called the high-precision A^A^ potentials. Those in coordinate 
space, some of which are shown in Fig. [T], reflect characteristic features of the A^A^ interaction at different 
lengthscalesfHElEllllE]: 

(i) The long range part of the nuclear force (at relative distance r > 2 fm) is dominated by one-pion 
exchange introduced by Yukawa[9j. Because of the pion's Nambu-Goldstone character, it couples 
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Figure 1: Three examples of the modern NN potential in ^5*0 (spin-singlet and S-wave) channel: CD- 
Bonn [6], Reid93[7] and Argonne UislH]. Taken from Ref. [T5] . 



to the spin-isospin density of the nucleon and hence leads to a strong spin-isospin dependent force, 
namely the tensor force. 

(ii) The medium range part (1 fm < r < 2 fm) receives significant contributions from the exchange 
of two-pions (tttt) and/or heavy mesons (p, u, and a). In particular, the spin-isospin independent 
attraction of about 50 - 100 MeV in this region plays an essential role for the binding of atomic 
nuclei. 

(iii) The short range part (r < 1 fm) is best described by a strong repulsive core as originally introduced 
by Jastrow [lOj. Such a short range repulsion is important for the stability of atomic nuclei 
against collapse, for determining the maximum mass of neutron stars, and for igniting the Type 
II supernova explosions [HI [12], [13] . 

It is then a challenge for the theoretical particle and nuclear physics communities to extract these 
properties of the A^A^ interaction from first principle non-perturbative QCD calculations, in particular 
lattice QCD simulations. A theoretical framework suitable for such a purpose was first proposed by 
Liischer[T3]. For two hadrons in a finite box with the size L x L x L with periodic boundary conditions, 
an exact relation between the energy spectra in the box and the elastic scattering phase shift at these 
energies was derived. If the range of the hadron interaction R is sufficiently smaller than the size of 
the box R < L/2, the behavior of the two-particle Nambu-Bethe-Salpeter (NBS) wave function ipi^) 
in the interval -R < |r| < L/2 has sufficient information to relate the phase shift and the two-particle 
spectrum. Liischer's method bypasses the difficulty to treat the real-time scattering process on the 
Euclidean lattice. Furthermore, it utilizes the finiteness of the lattice box effectively to extract the 
information of the on-shell scattering matrix and the phase shift. 

A closely related but an alternative approach to the A^A^ interactions from lattice QCD has been 
proposed recently [IHl |16l [I7j. The starting point is the same NBS wave function V^(r) as discussed 
in Ref. [H]. Instead of looking at the wave function outside the range of the interaction, the authors 
consider the internal region \r\ < R and define an energy-independent non-local potential U{r, r') from 
'?/'(r) so that it obeys the Schrodinger type equation in a finite box. Since U{r, r') for strong interactions 



is localized in its spatial coordinates due to confinement of quarks and gluons, the potential receives 
only weak finite volume effect in a large box. Therefore, once U is determined and is appropriately 
extrapolated to L — )■ cxd, one may simply use the Schrodinger equation in infinite space to calculate 
the scattering phase shifts and bound state spectra to compare with experimental data. A further 
advantage of utilizing the potential is that it is a smooth function of the quark masses so that it is 
relatively easy to handle on the lattice. This is in sharp contrast to the scattering length which shows 
a singular behavior in the quark mass corresponding to the formation of the NN bound state. 

Since the recent progress for the study of the NN interaction by the first method has already been 
reviewed in Ref. |18j. the recent progress for the second method is mainly considered in this review. In 
Sec. [21 the strategy of Ref. [151 [HI HZI to define the NN potential in QCD is explained in detail, and the 
lattice formulation is introduced in Sec. [31 Results of lattice QCD calculations for A^A^ potentials are 
given in both quenched and full QCD in Sec. [H Central potentials at the leading order of the velocity 
expansion is shown to reproduce qualitative features of the A^A^ potential such as the repulsion at short 
distance and the attraction at medium to long distances. The tensor potential, which exists also at 
leading order, is extracted. Contrary to the case of the central potentials, it does not have a repulsive 
core. Higher order contributions in the velocity expansion are also investigated and shown to be small 
at low energy and low orbital angular momentum L. In Sec. [SI the method to extract the potential is 
applied to the hyperon-nucleon interactions such as NS and A^A systems. Interactions between octet 
baryons in general are also investigated in the flavor SU(3) limit, where up, down and strange quark 
masses are all equal. In Sec. [6l we also consider a recent attempt to understand the origin of the 
repulsive core in the A^A^ potential. Using the operator product expansion and renormalization group 
analysis in QCD, the potential derived from the NBS wave function in Sec. [His shown to have a repulsive 
core, whose functional form is also theoretically predicted. In Sec. U\ two extensions of the potential 
method are considered, together with explicit applications of these extensions to hadron interactions. 
One is the extension of the potential method to inelastic scattering, in order to investigate the AA 
system, while the other is the extraction of the potential from the time dependent NBS wave function 
in lattice QCD. With the latter method, the existence of the H-dihaijon is investigated in the flavor 
SU(3) limit. In Sec. [H applications of the method to the three nucleon force, meson-baryon potentials 
and the potential in 2-color QCD are considered. Brief concluding remarks are given in Sec. [91 

2 Strategy to extract potential in QCD 

2.1 Nambu-Bethe-Salpeter (NBS) wave function and its asymptotic be- 
havior 

A key quantity to extract the potential from QCD is the equal time Nambu-Bethe-Salpeter wave 
function, defined by 

y.^(x)e-^* = {0\T{N{r + ^,t)Nir,t)}\2N,W,s,S2), (1) 



where \2N,W, S1S2) is an eigen-state of QCD for two nucleons with total energy W = 2J'k'^ + m% 
and the total three-momentum P = 0, whose helicities are denoted by si, S2- Here the local nucleon 
operator is given by 

Na{x) ^ (^^^|^])=e'^''^K(x)C754(x))g«(x), 9(^) = ( JJ(tJ ) , (2) 

where x = (x, t), the color indices are denoted by a,b,c, and a is the spinor index. The charge 
conjugation matrix in the spinor space is given by C = 7274, and p, n denotes proton and neutron 



operators while u, d denote up and down quark operators. Note that (p^ imphcitly has two pairs of 
spinor-flavor indices which come from Na{v + yi,t)Np{Y,t) and two hehcity indices Si and S2- 

The most remarkable property of the above NBS wave function is explained as follows. If the W is 
smaller than the threshold energy for one-pion production {i.e. W < 2mN + ttl.^ ), then its asymptotic 
behavior for large |x| can be evaluated |19[ [T7] . The helicity component in the spin singlet channel 

(5* = 0) is given by |siS2) = "75 ( I "I" I' ~'~^) + I ~ I' ~|))' "where the relative + sign is our convention. 
For this case, we have 

^»'(r),.„ == Z Z'-''(S = 0)y,,(n,) ""<^'- - ";f + ^'°"^)) e-'..<» (3) 

where r = |r|, fc = |k|, and 6is{k) is the A^A^ scattering phase shift in QCD with the total angular 
momentum I and the total spin S, which is determined by the unitarity of S-matrix in QCD below the 
inelastic threshold[17]. Here Y/m(^r) is the spherical harmonic function with the solid angle Qr of r- 
The coefficient Z''^(0) has spinor components a,/3 and is given by 

Zi'^(O) = ZDUn^)U^^i\/)U^^i-V)x^^iO,0) (4) 

where Z is the wave function renormalization for the nucleon operator N{x), Dl^^ is the Wigner D- 
matrix, and Uaai"^) ^^^ t/go(— V) are the 4x2 matrices acting on the 2x2 matrix Xasi^y ^z)- Explicitly 
we have 

-ia -V 



U(V) = V»' + "«(,/.x., j^T^^j (5) 

X(0,0) = -^i<T2, X(1,0) = ;^<T1, X(l,±l) = i(/2X2±<T3). (6) 

For the spin triplet channel {S = 1), the asymptotic behavior of (p^ is more involved but schematically 
written as 

^•"(r),., oci:^i..(».) ""'*'"';f'''""" e-'..l-). (7) 

An explicit form of the asymptotic behavior for the NBS wave function in the triplet channel is given 
in appendix |X1 

The asymptotic behaviors in eqs. ([3]) and ([7]) tell us that the NBS wave function at large separation 
r describes the scattering wave of the quantum mechanics whose phase shift agrees with the phase of 
the S-matrix in QCD. Therefore the NBS wave function satisfies the free Schrodinger equation at large 
r as 



Ho 

2/x " 



^^(r) ^ 0, Ho ^' 



2/i 

at W < 2m N + rriT,, where /i = 171^/2 is the reduced mass of the A^A^ system. Note that these 
properties hold without using the non-relativistic approximation/expansion. In particular, only the 
upper components of the spinor indices for the NBS wave function (a = 1, 2 and /3 = 1, 2) are enough 
to reproduce all A^A^ scattering phase shifts Sis{k) with / = 0,1,2,3, ■■■ and S = 0,1. From these 
properties, (the upper spinor components of) the NBS wave function can be regarded as the "wave 
function" of the A^A^ system a.t W < 2mN + rn.,^. 

It is also noted that the equal-time constraint for the NBS wave function here is not a restriction to 
the extraction of physical observables such as the scattering phase shift, as evident from the fact that all 
informations on the scattering phase shifts are encoded in the asymptotic behavior of the equal-time NBS 
wave function. Moreover, the equal-time NBS wave function with non-zero total momentum (P ^ 0) 
is equivalent to the NBS wave function in the center of mass frame with non-zero time separation. In 
this case the 4-dimensional distance between two nucleon operators is always space-like. 



2.2 Non-local potential from the NBS wave function 

Since the NBS wave function satisfies the free Schrodinger equation at large r, one can define short- 
ranged non-local potential by 



It is noted that the spinor indices a, /3, 7, 6 here run from 1 to 2, since all A^A^ scattering phase shifts can 
be reproduced from the NBS wave function with a,/? G {1,2} as discussed in the previous subsection. 
Therefore Uap-^^s has 4x4 components, which can be determined from 4 components of (p^o for 4 different 
combinations of (si,S2). Note that since the NBS wave function ip^ is multiplicatively renormalized, 
the potential f/(x, y) is finite and does not depend on the particular renormalization scheme . 

The non-local function f/(x, y) is shown to be energy-independent as follows. Let Vf^ be the space 
spanned by the wave function with W < Wth = 2mjv + rriT,: V^h = {v^'^lVT < PVth}, and the projection 
operator to V^h is defined as 



JWt 9<Wth 

p{W^)dW^P{Wu^,y) (10) 



Wi<W,i, 



where p{W) is the density of states at energy W, and A^ ^(Vri,PV2) is the inverse of the hermitian 
operator N{Wi, W2) defined by 



N{Wi, W2) = / ^^Hr)V^'(r) d\ (11) 

so that 

[ p{W)dW N{Wi,W)N~\W,W2) = -^6{Wi-W2). (12) 

J p{Wi) 

The non-local potential is then defined by 

^^'Hx,y) = / p{Wi)dWip{W2)dW2[Ek-Ho]v'^'i^)N~\W,,W2)v'^'{y)^ 

= f p{W,)dW,[Ek-Ho]P{W,;^,y). (13) 

It is easy to see that the above non-local potential satisfies eq.(l9]) at W^ < Wth as follows. 

/f/(x,y)^-¥,^(y)A = / piW^)dW^[E,-H,]^'^^i^)^6iW,-W) 

= ^(PVth - W) [Ek - Ho] v9^(x). (14) 

It should be noted that the non-local potential f/(x, y) which satisfied eq.([9]) at ly < Wth is not unique: 
For example, we can add the following term 

/ p{W)dWfw{^)P{W-^,y) (15) 



with arbitrary functions /vk(x) to the non-local potential [/(x, y) without affecting eq.Q at ly < PVth- 
The non-local potential f/(x, y) in eq. flT3|) is energy independent by construction. 
Alternatively we can define a different non-local potential by 

t/°°(x,y)=/ p{W)dW[Ek-Ho]P{W;^,y), (16) 

Jo 



which satisfies eq. iQ for all W. This potential, however, becomes long-ranged, due to the presence of 
inelastic contributions above Wth- The extension of this method to non-elastic cases will be discussed 
in Sec. [71 

In Ref. [18], it is claimed that the NBS wave function satisfies the Schrodinger equation with a 
non-local and energy-dependent potential. In general this is true but, as shown here, there is a scheme 
which makes the non-local potential energy-independent. This is sufficient for the strategy considered 
in this report. 

2.3 Velocity expansion of the non-local potential 

In principle, if one knows all NBS wave functions (p^ , the non-local potential U can be constructed 
according to eq. (TT^ . In practice, however, one can obtain only a few of them corresponding to the 
ground state as well as a few low lying excited states in lattice QCD simulations. Therefore, for practical 
applications, it is convenient to expand the non-local potential in terms of the velocity(derivative) with 
local functions as 

t/(x,y) = V(x,V)5=^(x-y). (17) 

At the lowest few orders we have 

V(r,V) = Voir) + K(r)ai ■ ^2 + VT{r)Su + VLs{r)L ■ S +0(V^), (18) 

LO NLO 

where r = |r|, cij is the Pauli-matrix acting on the spin index of the i-th nucleon, S = [cxi + ct2)/2 is the 
total spin, L = r x p is the angular momentum, and 

(r-ffi)(r-(T2) ^ ^ . ^ 

Oi2 = S (Ti ■ 0-2 (,19) 

is the tensor operator. Each coefficient function is further decomposed into its fiavor components as 

Vx{r) = VUr) + V],{r)n-r2, X = 0, a,T, LS, ■ ■ ■ , (20) 

where Tj is the Pauli-matrix acting on the fiavor index of the i-th nucleon. The form of the velocity 
expansion flTS]) agrees with the form determined by symmetries (20 j. 

For the leading order of the velocity expansion, the local potential is given by 

V^L«(r) = Voir) + K(r)ai ■ c?2 + VV(r)^i2, (21) 

which can be obtained from the NBS wave function at one value of W. Since S12 = for the spin single 
state, for example, we have 

Vir, S = 0)^ Vir) - 3K(r) = ^^^ ~ ^°^ f ^^^ (22) 

2.4 Remarks 

There are several remarks on the extraction of the potential from QCD described in the previous 
sections. 

First of all, it should be noted that the potential itself is not a physical observable, and it is therefore 
not unique. In particular the potential depends on the choice of the nucleon operator to define the NBS 
wave function. The local nucleon operator is one choice here but other definitions are equally possible, 
though the local operator is a convenient choice for the reduction formula of composite particles such as 

7 



the nucleon[2Tl [221 [23]. A choice of the nucleon operator to define the NBS wave function is considered 
to be a "scheme" to define the potential The potential is therefore a scheme dependent quantity, while 
physical observables such as the scattering phase shift and the binding energy of the deuteron are of 
course scheme independent. 

Is such a scheme-dependent quantity useful ? The answer to this question is probably "yes", since 
the potential is useful to "understand" or "describe" the phenomena. For example, the repulsive core 
best summarizes the behavior of the NN scattering phase shift at larger energy in terms of the short 
distance behavior of the potential. A well-known example for a scheme- dependent but useful quantity 
is the running coupling in quantum field theory. Although the running coupling is of course scheme 
dependent, it is useful to understand the deep inelastic scattering data at high energy {i.e. asymptotic 
freedom) . 

Among different schemes, there exist of course good schemes. A good convergence of the perturbative 
expansion for a certain class of observables is a reasonable criteria for a good scheme in the case of 
the running coupling. In the case of the potential, a good convergence of the velocity expansion is 
important. In this respect, a completely local and energy-independent potential is the best one, and 
moreover the inverse scattering method tells us that it must be unique if it exists. We think that such 
a local and energy-independent potential exists only if no inelastic scattering appears at all values of 
energy {i.e. Wth = oo). 

There exists a criticism that the repulsive core in the phenomenological potentials is meaningless 
since the low energy scattering data can be described equally well by other potentials without the 
repulsive core. In other words, the repulsive core can be removed by a unitary transformation of the 
wave function without changing the scattering phase shift at low energy. Even though this may be 
true, the resulting "potential" is highly non-local, and therefore is beyond the concept of a potential. 
This criticism thus corresponds to claiming that the asymptotic freedom is meaningless since it can be 
removed by some non-perturbative scheme of the "coupling", which is however beyond the concept of 
a coupling. 

For the definition of the potential in QCD from the NBS wave function, the convergence of the ve- 
locity expansion can be checked by examining the energy {W) dependence of the lower order potentials. 
For example, if we have cp^'^ for ra = 1, 2, ■ ■ ■ A^, we can determine the A^ — 1 unknown local functions of 
the velocity expansion in A^ different ways. The variation among A^ different determinations gives an 
estimate of the size of the higher order terms neglected. Furthermore one of these higher order terms 
can be determined from ip^" for n = 1,2,- ■■ N. The convergence of the velocity expansion will be 
investigated explicitly in Sec. [H 

Let us consider what we have shown so far from a slightly different point of view. Instead of adopting 
the point of view that the "potential" can be defined in QCD, the analysis in this section shows that 
the use of quantum mechanics with potentials to describe the NN scattering can be justified in a more 
fundamental quantum field theory (QCD) through the NBS wave function, whose asymptotic behavior 
encodes phases of the S-matrix for the NN scattering. If the local approximation for the potential 
defined from the NBS wave function is good in the low energy region, one can use the local potential 
combined with quantum mechanics to investigate nuclear physical! . 

It is now clear that there is no unique definition for the NN potential. Ref . [HI [251 126] , however, 
criticized that the NBS wave function is not "the correct wave function for two nucleons" and that its 



"'^It is more legitimated to show if three or more body potentials remain subdominant in the same framework in QCD 
before using the potential in nuclear physics. There exists another approach to investigate nuclei, by combining the 
chiral effective field theory with lattice calculations. (See Ref. [24] and references therein.) The chiral effective can easily 
incorporate may-body interactions order by order in the chiral expansion, though there are many free parameters to 
be determined from experimental data. On the other hand, the potential from lattice QCD does not contain such free 
parameters. Therefore it is interesting and important to determine some parameters of the chiral effective theory from 
the potential in lattice QCD. 



relation to the correct wave function is given by 

(^^(r) = Zjv^(|r|)(0|T{iVo(x + r, 0)iVo(x, 0)}|2iV, W, s,, S2) + ■ ■ ■ (23) 

where No{'x,t) is "a free- field nucleon operator" and the ellipses denote "additional contributions from 
the tower of states of the same global quantum numbers". Thus (0|T{A^o(x+r, 0)A^o(x, 0)}|2A^, W, si, S2) 
is considered to be "the correct wave function". In this claim it is not clear what is "a free-field nucleon 
operator" in the interacting quantum field theory. An asymptotic in or out field operator may be a 
candidate. If the asymptotic field is used for Nq, however, the potential defined from the wave function 
identically vanishes for all r by construction. To be more fundamental, a concept of "the correct wave 
function" is doubtful. If some wave function were "correct", the potential would be uniquely defined 
from it. This clearly contradicts the fact discussed above that the potential is not an observable and 
therefore is not unique. This argument shows that the criticism of Ref. [THJ |25l [26] is fiawed. 

3 Lattice formulation 

In this section, we discuss the procedure to extract the NBS wave function from lattice QCD simulations. 
For this purpose, let us consider the correlation function on the lattice defined by 

F{r,t-to) = (0|T{iV(x + r,t)iV(x,t)}:7(to)|0) (24) 

where J^ito) is a source operator which creates a two nucleon state of states and its explicit form will 
be considered later. By inserting a complete set and considering baryon number conservation, we have 

F(r,t-to) = (0|T{iV(x + r,t)iV(x,t)} ^ \2N,Wn,SuS2){2N,Wn,s,,S2\J{tom 

n,si,S2 

= J2 ^n,si,s.¥''^"(r)e-^"(*-*«), A„,,„,, = (2iV,iy„,Si,S2|J(0)|0). (25) 



n,s\, S2 



For a large time separation (t — tg) — ?■ 00, we have 

lim F(r,t - to) = Ao^^°(r)e-^°(*-*«) + 0(e-^"^«(*-*°)) (26) 

{t—to)—>oo 

where Wq is assumed to be the lowest energy of NN states. Since the source dependent term Aq is just 
a multiplicative constant to the NBS wave function if^°{r), the potential defined from ip^°{r) in the 
previous section is manifestly source-independent. Therefore the statement that the potential in this 
scheme is " source- dependent" in Ref. [27] is clearly wrong. 

In this extraction of the wave function, the ground state saturation for the correlation function F in 
eq. f l26p is important. In principle, one can achieve this by taking a large t — to. In practice, however, 
F becomes very noisy at large t — to, so that the extraction of ip^° becomes difficult at large t — to- 
Therefore it is crucial to find the region of t where the ground state saturation is approximately satisfied 
while the signal is still reasonably good. The choice of the source operator becomes important in order 
to have such a good t-region. 

3.1 Choice of source operator 

We can choose the source operator jf to fix quantum numbers of the state |2A^, W, Si, S2) such as {J, J^)- 
Since a lattice QCD simulation is usually performed on the (finite) hyper-cubic lattice, we consider the 
cubic transformation group S0{3, Z) instead of the S0{3, R) as the symmetry of 3-dimensional space. 
Therefore the quantum number is classified in terms of the irreducible representation of S0{3,Z), 

9 



Table 1: The number of each representation of S0{3,Z) which appears in the angular momentum L 
representation of SO {3, R). P = (— 1)^ represents an eigenvalue under parity transformation. 



L P 


Ai 


A2 


E Ti 


T2 


0(S) + 


1 











1(P) - 








1 





2(D) + 








1 


1 


3(F) - 





1 


1 


1 


4(G) + 


1 





1 1 


1 


5(H) - 








1 2 


1 


6(1) + 


1 


1 


1 1 


2 



Table 2: The decomposition for a product of two irreducible representations, Ri 
representations in 5*0(3, Z). Note that Ri® R2 = R'i® Ri by definition. 



R2, into irreducible 





Ai 


A2 


E Ti T2 


Ai 


Ai 


A2 


E Ti T2 


A2 


A2 


Ai 


E T2 Ti 


E 


E 


E 


Ai®A2®E Ti © T2 Ti © T2 


Ti 


Ti 


T2 


Ti © r2 Ai © E © Ti © T2 A2®E®Ti®T2 


T2 


T2 


Ti 


Ti © T2 yl2 © E © Ti © T2 Ai © E © Ti © T2 



denoted by Ai, A2, E, Ti, T2 whose dimensions are 1, 1, 2, 3, 3. A relation of irreducible representations 
between SO {3, Z) and SO {3, R) is given in tabled] for L < 6, where L represents the angular momentum 
for the irreducible representation in 5*0(3, R). For example, the source operator jf(to) in the Ai 
representation with positive parity generates states with L = 0, 4, 6, ■ ■ ■ at t = to, while the operator in 
the Ti representation with negative parity produces states with L = 1, 3, 5, ■ ■ ■. For two nucleons, the 
total spin 5 becomes l/2®l/2 = 1©0, which corresponds to Ti(5 = 1) and Ai{S = 0) of the 50(3, Z). 
Therefore, the total representation J for a two nucleon system is determined by the product -Ri © -R2, 
where Ri = Ai,A2,E,Ti,T2 for the orbital "angular momentum" while R2 = Ai,Ti for the total spin. 
In table |21 the product -Ri © R2 is decomposed into the direct sum of irreducible representations. 
Most of the results in this report use the wall source at t = to defined by 



J-^'%)a,J, = NZI%)N^^^%) 



(27) 



where a, /3 = 1,2 are upper component spinor indices while f,g are flavor indices. Here A^™'^^^(to) is 
obtained by replacing the local quark field q{x) of N{x) by the wall quark field. 



q-^\to) = Y.q{^,to 



(2J 



with the Coulomb gauge fixing only at t = to- Note that this gauge- dependence of the source operator 
disappears for the potential. All states created by the wall source have zero total momentum. Among 
them the state with zero relative momentum has the largest magnitude. The most important reason to 
employ the wall source is that the ground state saturation for the potential at long distance is better 
achieved for the wall source than other sources such as the smeared source. 

By construction, the source operator jf™'^^^(to) has zero orbital angular momentum at t = to, which 
corresponds to the Ai representation with positive parity. By the spin projection operator P^^\ e.g. 
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piS-0) — (j^ QT^(^ p{S-i,s^-o) — u_^^ ^Q gx the J of the source as 

J{to; J^=+, /) = Pl^'Jj-^\toU,fg (29) 

where P = ± is the parity and / = 1, is the total isospin of the system. Since the nucleon is a fermion, 
an exchange of two nucleon operators in the source should give a minus sign. This fact fixes the total 
isospin once the total spin is given: {S, I) = (0, 1) or (1, 0). (Note that S,I = are antisymmetric while 

5.1 = 1 are symmetric under the exchange.) Since Af ^ Ai{S = 0) = Af and Af (g) Ti{S = 1) = T{^, 
the state with either {J^,I) = (Af,!) for the spin-singlet or {J^,I) = (T{^,0) for the spin-triplet is 
created at t = to by the corresponding source operator. The NBS wave function extracted at t > to 
has the same quantum numbers ( J^, /) as they are conserved under QCD interactions. In addition the 
total spin 5* is conserved at t > to for the two-nucleon system with equal up and down quark masses: 
Under the exchange of the two particles, the constraint (— l)'5'+i+^+ip = _ i must be satisfied due to 
the fermionic nature of nucleon, while the parity P and the isospin / are conserved in this system. 
Therefore 5* is conserved. It is noted, however, that L is not conserved in general. While the state with 
( J^, /) = {Af, 1) always has L = Af even at t > to, the one with ( J^, J) = (Tf', 0) has both L = Af 
and L = E^ componentcl at t > to, which corresponds to L = and L = 2 in 5*0(3, R), respectively. 
Note that J or L in this report is used to represent the total or orbital quantum number of 50(3, Z) 
as well as 5*0(3, R), depending on the context. 

The orbital angular momentum L of the NBS wave function can be fixed to a particular value by 
the projection operator P^^' as 

¥.^(r; J^, /, L, S) = p(^)p(^)y,^(r; J^, /) (30) 

where (p^{r] J^, I) is extracted from 

P(r, t - to; J^, /) =. A( J^, /)^^(r; J^, J)e-^(*-*«), A( J^, /) = (2iV, W\J{0; J^, /) |0) (31) 

for large t — to. Here we also apply the total spin projection operator P^^^ but this is redundant since 
the total spin 5, already fixed by the source, is conserved as mentioned above. The projection operator 
P^^'^ to an arbitrary function 0(r) is defined in general by 

P^'W)^^ E x'{9)<l^{9-'-r) (32) 

geSO{3,Z) 

for L = Ai,A2,E,Ti,T2, where x^ denotes the character for the representation L of the cubic group 
5*0(3, Z), g is one of 24 elements in 5*0(3, Z) and di is the dimension of L. 

3.2 Leading order potential: spin-singlet case 

We present the procedure to determine potentials at the reading order (LO): 

VLO(r) = Voir) + K(r)(ai • 0^2) + VT{r)Su. (33) 

Since 5*i2 = and a\ ■ 02 = —3 for the spin-singlet case, the LO central potential in the case of 
[J^, I) = {Af, 1) state is extracted as 

T/ f^^(5,/)=(o,i) ^ v'=^(r) - 3V'='(r) - ^^' " ^°^ '^'^^''' ^^' ^ = '^,^1,3 = 0) . . 

Vc[r) _l/o [r) 31/, [r) - ^w^^. ^f,! = 1,A,, S = 0) ' ^^' 



^This can be seen from Table [5] for R2 — Ti (spin-triplet), which also tells us existences of i = Tf and L = Tf 
components in addition. These extra components are expected to be small since they appear as a consequence of the 
violation of SO{3, R) on the hyper-cubic lattice. 
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where V^^^ = V^ + VJ in isospin space. The potential Vc{'r)^^'^''^^^'^^ in the above is often referred 
to as the central potential for the ^5*0 state, where the notation '^^'^^Lj represents the orbital angular 
momentum L (see table [T]), the total spin 5* and the total angular momentum J of J = L + S. It is 
noted, however, that in the leading order of the velocity expansion, the potential does not depend on the 
quantum number of the state J = L = Ai. Moreover the Ai state may contain L = 4, 6, ■ ■ ■ components 
other than L = 0, though the L = component may dominate in the state. Therefore it is more precise 
to refer to Vdn^)^^'^^^^^'^'' as the spin-singlet (isospin-triplet) central potential determined from the state 
with J = L = Ai. A possible difference of spin-singlet central potentials between this determination 
and others such as the one determined from J = L = E gives an estimate for contributions from higher 
order terms in the velocity expansion. 

3.3 Leading order potential: spin-triplet case 

Both the tensor potential Vr and central potential Vc appear at the LO for the spin-triplet case. Let 
us consider the determination from the {J^, I) = {T^, 0) state. The Schrodinger equation for this state 
becomes 

Ho + l^c(r)(^'^)=(^'°) + VT{r)S^,] ^^(r; J^ = T+, J = 0) = E,^^(r; J^ = T+, / = 0) (35) 

where the spin-triplet central potential is given by 

V^^(,)(^./)=(i.o) ^ vt%r) + V^'ir), V^' = V^ - 3V^. (36) 

The projections to Ai and E components read 

V^Z = P(^^Vr,(r;J^ = T+,/ = 0) (37) 

Q^Z ^ P(^Vr,(r;J^ = T+,/ = 0)^(l-P(^^))^r,(r;J^ = T+,/ = 0). (38) 

The last quantity in eq. (!38|) is an approximation of the first line and a difference comes from Ti and T2 
components, which are expected to be small. This approximate representation for Q is often employed 
in numerical simulations. 

Using these projections, Vc and Vr can be extracted as 

VcirY'''^ = i5fc-^([Q5i2^^]a/.(r)ifo[^^^]a/3(r)-[P5i2^^]„/3(r)i/o[Q^^]./3(r)) (39) 

Vrir) = ^ {[Q^'^U{r)Ho[V^^U{r) - [V^^U{r)Ho[Q^^U{r)) (40) 

A(r) ^ [QSuv'^Uir)[P^'^Uir)-[PSu^^Uir)[Q^'^Uir). (41) 

In numerical simulations, (a,/3) = (2, 1) is mainly employed. 

If one neglects Vr by putting V^ = in the above, one obtains the effective central potential for the 
spin-triplet (isospin-singlet) state as 

The difference between Vc and Vc^ is 0(y^) in the second order perturbation for small Vt- 
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3.4 A comparison with the finite volume method in lattice QCD 

In this subsection we briefly compare the potential method with the direct extraction of the phase shift 
via the finite volume method in lattice QCD. 

First of all, by construction, the potential method gives the correct phase shift at k = JW^ /A — mj^ 
where W is the total energy of the state from which the NBS wave function is defined, while phase 
shifts at other values of k are approximated ones obtained in the velocity expansion of the non-local 
potential. 

Secondly, the finite size correction to the potential is expected to be small. Indeed the finite volume 
method for the extraction of the phase shift in lattice QCD assumes that there is no finite size correction 
to the potential as long as the volume is large enough so that the interaction range of the potential is 
smaller than half of the lattice extension, L/2. Under this condition, there exists an asymptotic region 
in the periodic box where the scattering wave satisfies the free Schrodinger equation with a specific 
value of the energy, from which one can determine the phase shift at certain values of k in the infinite 
volume. This is Liischer's finite volume formula for the phase shift (Hj . 

Thirdly, we also expect that the quark mass dependence of the potential is much milder than that 
of physical observables such as the scattering length. While the scattering length is almost zero at the 
heavy quark mass region, it diverges when the bound state is formed at the lighter quark mass region. 
In this situation, the scattering length varies from zero to infinity as the quark mass changes [28]. Such 
a drastic change of the scattering length can easily be realized by a small change to the shape of the 
potential as a function of the quark mass. 

Let us assume that higher order terms in the velocity expansion give negligible contributions at 
low energy so that the leading order local potential well reproduces the scattering phase shift. In this 
situation, some problems of the finite size method can be avoided by using the potential method. To 
extract the phase shift in the finite size method in lattice QCD, one has to assume that one particular 
angular momentum gives a dominant contribution among possible angular momenta in a given rep- 
resentation of the cubic group. For example, although a state in the Ai representation contains not 
only an L = contribution but also L = 4, 6, ■ ■ ■ contributions, one usually assumes that the L = 
contribution dominates so that the energy shift in the finite volume is related to the scattering phase 
for the L = state. In the case of the potential, on the other hand, such an assumption is unnecessary. 
One can determine the local potential in the velocity expansion from the Ai state without specifying the 
dominant angular momentum. Once the potential is obtained, one can calculate the scattering phase 
shift for an arbitrary L by solving the Schrodinger equation in the infinite volume with the extracted 
potential. Furthermore, one can check the assumption made for the finite size method by comparing 
sizes of the scattering phases among different L's. 

4 Results for nuclear potentials from lattice QCD 



4.1 Quenched QCD results for (effective) central potentials 

Let us show results in the quenched QCD, where creations and annihilations of virtual quark-antiquark 
pairs are all neglected. For the simulations, the standard plaquette gauge action is employed on a 
32"^ lattice at the bare gauge coupling constant /3 = Q/g"^ = 5.7, which corresponds to the lattice 
spacing a ^ 0.137 fni (1/a = 1.44(2) GeV), determined from the p meson mass in the chiral limit, 
and the physical size of the lattice L ~ 4.4 fm[T5]. As for the quark action, the standard Wilson 
fermion action is employed at three different values of the quark mass corresponding to the pion mass 
m^ ~ 731, 529, 380 MeV and the nucleon mass tun — 1560, 1330, 1200 MeV, respectively. 
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Figure 2: (Left)The NN wave function for the spin-singlet and spin-triplet channels in the orbital Af 
representation at m^ ~ 529 MeV and a ~ 0.137 fm in quenched QCD. The insert is a three-dimensional 
plot of the spin-singlet wave function ip^{x,y,z = 0). (Right) The NN (effective) central potential for 
the spin-singlet (spin-triplet) channel determined from the orbital Af wave function. Both figures are 
taken from Ref. [17] . 



Fig. [2]^Left) shows the NBS wave functions for the spin-singlet and the spin-triplet channels in the 
orbital Ai representation at m^ ~ 529 MeV. These wave functions are normalized to be 1 at the largest 
spatial point r ~ 2.2 fm. 

The central potential in the spin-singlet channel and the effective central potential in the spin- 
triplet channel reconstructed from the wave functions at m.,^ ~ 529 MeV are shown in Fig. [2](Right). 
These potentials reproduce the qualitative features of the phenomenological NN potentials, namely the 
repulsive core at short distance surrounded by the attractive well at medium and long distances. From 
this figure one observes that the interaction range of the potential is smaller than 1.5 fm. Therefore the 
box size L ~ 4.4 fm is large enough to extract the phase shift by the finite size method, and furthermore 
the finite size corrections to the potentials themselves are expected to be small. Labels ^5*0 and ^Si of 
the potentials in the figure represent the fact that potentials are determined from Ai wave functions, 
which are dominated by S wave components. 

Instead of calculating the energy shift due to the finite size, one can extract the asymptotic momen- 
tum k, by fitting the NBS wave function if{r) at large distance with the Green's function G{r; k"^) in a 
finite and periodic box for the Helmholtz equation (V^ -|- k'^)G{r; k"^) = —6iat{r) with 6iat{r) being the 
periodic delta-function. Explicitly it is given by 



-^ gi(27r/L)n-r 

G(r;P) = ^E^^2vr/L)%^-A; 



2' 



(43) 



The asymptotic momentum k is related to the scattering phase shift 5o(^) or the scattering length oq 
for the S-stateqfl as 

kcot6o{k) = -i-Zoo(l;g') = - + 0{e), (44) 

a/ttL ao 

where Zqq{1, q^) with g = 1^ is (the analytic continuation of) the generalized zeta-function Zqq{s, q^) = 
y^ (n^ — q^)^^ ■ Fig. [3] shows the fits of the wave functions in the interval 11a ~ 1.5 fm < r < 

neZ3 



We here assume that the dominant component of the scattering wave in Ai representation has L = 0. 
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Figure 3: (Left) The fit of the A^A^ wave functions at 171^^ ~ 529 MeV for the spin-singlet channel in the 
orbital Af representation using the Green's function in the fit range 11 < r/a < 15. (Right) A similar 
fit for the spin-triplet channel. Taken from Ref. [T7]. 



Table 3: Effective center of mass energies E obtained from the asymptotic momenta and the scattering 
length oq at different pion masses. Taken from Ref. [17J. 





E[MeV] 


flo [fm] 


m^[MeV] 


spin-singlet spin-triplet 


5*0 5*1 


731.1(4) 
529.0(4) 
379.7(9) 


-0.40(8) -0.48(10) 
-0.51(9) -0.56(11) 
-0.68(26) -0.97(37) 


0.12(3) 0.14(3) 
0.13(3) 0.14(3) 
0.15(7) 0.23(10) 



16a ~ 2.2 fm using the above form at m.^ = 529 MeV. This leads to the values of the effective energy 
E = k'^/rriM, which can be translated to the scattering length oq by the Liischer's formula (Hll). 

In FigJH we compare the A^A^ central potentials in the spin-singlet channel for three different pion 
masses. As the pion mass decreases, the repulsion at short distance and the attraction at medium 
distance are enhanced simultaneously. In table [31 we give values of E and the S-wave scattering length 
Oq, which show a net attraction of the NN interactions in both channels at these pion masses, though 
the absolute magnitudes of the scattering length ao are much smaller than the experimental values at 
the physical pion mass m^ ~ 140 MeV: Oq ^ = {'Sq) ~ 20 fm and Oq ^ = (^"S*!) ~ —5 fm. 

The above discrepancy is partly caused by the heavier pion masses and the absence of the dynamical 
quarks in quenched simulations. If we get closer to the physical pion mass in full QCD simulations, 
there should appear the "unitary region" where the A^A^ scattering length shows the singularity asso- 
ciated with the formation of the di-nucleon bound state, so that oq changes sign[28]. Therefore the 
A^A^ scattering length must become a non-linear function of the pion mass in this region. Unlike the 
scattering length, on the other hand, the A^A^ potential does not necessarily have singular behavior 
in the unitary region, as demonstrated in the well-known quantum mechanical examples such as the 
low-energy scattering between ultracold atoms. To check this in QCD, it is of course important to study 
the A^A^ potential in the full QCD simulations at lighter pion masses. 

In addition to the above reasoning, there is a possibility that extracted values of k'^ have large 
systematic uncertainties caused by the contamination of the excited states at large distance for the 
wave functions. 
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Figure 4: The central potentials for the spin-singlet channel from the orbital Af representation at three 
different pion masses in quenched QCD. Taken from Ref. jl7] . 
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Figure 5: (Left) (a, (3) = (2, 1) components of the orbital Af and non-A^i^ wave functions from J^ = T^ 
(and Jz = Sz = 0) states at m.^ ^ 529 MeV. (Right) The same wave functions but the spherical 
harmonics components are removed from the non-A^*" part. Taken from Ref. |17] . 



These A^A^ scattering lengths extracted from the NBS wave function agree in sign but are much 
smaller in magnitude than the previous quenched results from the finite size method in smaller volume [29] . 
while they disagree even in sign with the recent full QCD results form the finite size method (See Ref. [18] 
and references therein.). 



4.2 Tensor potential 

In Fig. [St^Left), we show the Ai and non-Ai components of the NBS wave function obtained from the 
J^ = Ti (and Jz = Sz = 0) states at m^r — 529 MeV, according to eqs. (|37|) and f p8|) . The Ai wave 
function is multivalued as a function of r due to its angular dependence. For example, the (a, /3) = (2, 1) 
component of the L = 2 part of the non-Ai wave function is proportional to the spherical harmonics 
^2o(^,0) oc 3cos^6' — 1. Fig. |5]^Right) shows the non-Ai component divided by Y2o{9,(j)). It is clear 
that the multivaluedness is mostly removed, showing that the non-Ai component is dominated by the 
D {L = 2) state. 
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Figure 6: (Left) The central potential Vc{rY^'^' and the tensor potential Vt(^) obtained from the 
J^ = Ti NBS wave function, together with the effective central potential V^^ 



f (r)(i'°), at m^ ~ 529 

MeV. (Right) Pion mass dependence of the tensor potential. The lines are the four-parameter fit using 
one-pion-exchange + one-rho-exchange with Gaussian form factor. Taken from Ref. |17] . 



Shown in Fig. [6] (Left) are the central potential Vc'(r)*^^'*^^ and tensor potential Vt(?") together with 



the effective central potential V^ 



eff/ 



Ui,o) 



at the leading order of the velocity expansion as given in eqs. 



dMD, (HOD and (112]), respectively. 

Note that V^^{r) contains the effect of VTir) implicitly as higher order effects through processes 
such as ^5*1 — > ^Di -^ ^Si. At the physical pion mass, V§^{r) is expected to obtain sufficient attraction 
from the tensor potential, which causes the appearance of a bound deuteron in the spin-triplet (and 
fiavor-singlet) channel while an absence of the bound dineutron in the spin-singlet (and flavor-triplet) 
channel. The difference between Vc{r)'^^'^'> and V^^{r) in Fig. [6] (Left) is still small in this quenched 
simulation due to relatively large pion mass. This is also consistent with the small scattering length in 
the previous subsection. 

The tensor potential in Fig. [6] (Left) is negative for the whole range of r within statistical errors and 
has a minimum around 0.4 fm. If the tensor potential receives a significant contribution from one-pion 
exchange as expected from the meson theory, VT{r) is rather sensitive to the change of the pion mass. 
As shown in Fig. [6] (Right), it is indeed the case: Attraction of Vr(r) is substantially enhanced as the 
pion mass decreases. 

The central and tensor potentials obtained from lattice QCD are given at discrete data points. For 
practical applications to nuclear physics, however, it is more convenient to parameterize the lattice 
results by known functions. Such a fit for Vt(^) is given by the form of one-pion-exchange -|- one-rho- 
exchange with Gaussian form factors as 



VAr) 



hiil-e-'^'"f\l + 



+ 



mpT 



[rupr) 



+ &3(l-e- 



-647-2 \ 2 



1 + 



+ 



mj^r {rriT^rY ) r 



(45) 



where 61,2,3,4 are the fitting parameters while mir {irip) is taken to be the pion mass (the rho meson 
mass) calculated at each pion mass. The fit line for each pion mass is drawn in Fig. [6] (Right). It may 
be worth mentioning that the pion-nucleon coupling constant extracted from the parameter 63 in the 
case of the lightest pion mass (m^r = 380 MeV) gives (7^^/(47r) = 12.1(2.7), which is encouragingly close 
to the empirical value. 
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Figure 7: (Left) The spin-singlet central potential Vc'(^)''°'^'' obtained from the orbital Af channel at 
i? ~ 45 MeV (red solid circles) and at -E ~ MeV (blue open circles) in quenched QCD at m^ ~ 529 
MeV. (Right) The same potentials at -E ~ 45 MeV, obtained from the orbital Af representation (red 
open circles) and from the T^ representation (cray solid circles). Taken from Ref. [32] . 



4.3 Convergence of the velocity expansion 

The potentials are derived so far at the leading order of the velocity expansion. It is therefore important 
to investigate the convergence of the velocity expansion: How good is the leading order approximation 
? How small are higher order contributions ? If the non-locality of the NN potentials were absent, the 
leading order approximation for the potentials would give exact results at all energies. The non-locality 
of the potentials therefore becomes manifest in the energy dependence of the potentials. 

So far the LO potentials are extracted with periodic boundary conditions in the spatial directions 
for quark fields. This leads to the lowest values of the effective center of mass energy E almost zero. 
To study the energy dependence, the leading order local potentials at i? ~ 45 MeV, realized by anti- 
periodic boundary conditions in the spatial directions, are calculated in quenched QCD at m^ ~ 529 
MeV and L ~ 4.4 fm[30l [3T| [32] . In this case, 4 types of "momentum- wall" sources, defined by 



„wall 



(to) = E'?(^'^o)/(x) 



(46) 



are employed, where /(x) = cos((±a; ± y + z)n/L). Note that /(x) = 1 corresponds to the wall source 
used in the periodic boundary condition. These momentum-wall sources induce L = T2 as well as 
L = Af states. 

In Fig. [Tl^Left), the spin-singlet potential Vc{r)^^'^^^^'^'^^ obtained from the L = Ai state at i5 ~ 45 
MeV (red circles) is compared with that at i? ~ MeV (blue circles), while a comparison is made for 
the spin-triplet potentials in Fig. [HI Vc'(r)'^'^'^)"'^^'°)(left) and Vrir) (right). Good agreements between 
results at two energies seen in these figures indicate that higher order contributions are rather small in 
this energy interval. In other words, these local potentials obtained at -E ^ MeV can be safely used to 
describe the NN scattering phase shift in both spin-singlet and -triplet channels between E = MeV 
and ii^ = 45 MeV at this pion mass in quenched QCD. 

Non-locality of the potential may become manifest also in its angular momentum dependence, since 
the orbital angular momentum L = r x p contains one derivative. In Fig. [7| (Right), the spin-singlet 
potential Vc{rY^'^^^^'^'^'^ obtained from the L = T2 state, whose main component has L = 2, is compared 
to the one from the L = Af state, whose main component has L = 0. In this case local potentials are 
determined at the same energy, i? ~ 45 MeV, but different orbital angular momentum. Although the 
statistical errors are rather large in the case of L = T2, a good agreement between the two is again 
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Figure 9: (Left) 2+1 flavor QCD results for the central potential and tensor potentials at rriT^ 
MeV. (Right) Quenched results for the same potentials at 711^^ ~ 731 MeV. Taken from Ref. [35] 
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observed, suggesting that the L dependence of the potential is small for the spin-singlet state. 

By these comparisons, it is observed that both energy and orbital angular momentum dependences 
for local potentials are very weak within statistical errors. We therefore conclude that contributions 
from higher order terms in the velocity expansion are small and that the LO local potentials in the 
expansion obtained at ii^ ^ MeV and L = are good approximations for the non-local potentials at 
least up to the energy ii^ ~ 45 MeV and orbital angular momentum L = 2. 

Hereafter "potential" in this report means the local potential at the leading order, unless otherwise 
stated. 

4.4 Full QCD results 

Needless to say, it is important to repeat calculations of NN potentials in full QCD on larger volumes 
at lighter pion masses. The PACS-CS collaboration is performing 2 + 1 flavor QD simulations, which 
cover the physical pion mass [331 [31] • Gauge conflgurations are generated with the Iwasaki gauge action 
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Figure 10: (Left) 2+1 flavor QCD results for the spin-singlet central potential from the orbital Af 
channel at three values of the pion mass. (Right) Scattering phase shifts in ^5*0 channel from the 
corresponding lattice potential given in (Left), together with the empirical one. Taken from Ref. 



and non-perturbatively 0(a)-improved Wilson quark action on a 32'^ x 64 lattice. The lattice spacing a 
is determined from m.,^, tjik and uiq as a ^ 0.091 fm, leading to L ~ 2.9 fm. Three ensembles of gauge 
configurations are used to calculate A^A^ potentials at {mT,,m]\[) ~(701 MeV, 1583 MeV), (570 MeV, 
1412 MeV) and (411 MeV,1215 MeV )|35] . 

Fig. inU^Left) shows the A^A^ local potentials obtained from the PACS-CS configurations at E' ~ 
and m^r = 701 MeV, which is compared with the previous quenched results at comparable pion mass 
m^r — 731 MeV but at a ~ 0.137 fm, given in Fig. [9]J^Right). Both the repulsive core at short distance 
and the tensor potential become significantly enhanced in full QCD. The attraction at medium distances 
tends to be shifted to the outer region, while its magnitude remains almost unchanged. These differences 
may be caused by dynamical quark effects. For a definite conclusion on this point, however, a more 
controlled comparison at the similar lattice spacing is needed. 

In Fig. [TOT Left). the spin-singlet central potential Vc{r)^^'^^ determined from the orbital Ai channel 
is plotted at three pion masses, while the spin-triplet central potential Vc{rY^'^^ and the tensor potential 
Vrir) from the orbital Af — T2 couple channel are given in Fig. [TTJ As in the quenched QCD, the 
repulsive cores at short distance, the attractive pocket at medium distance and the strength of the 
tensor potential are all enhanced as pion mass decreases. 

The phase shifts of the A^A^ scattering for ^Sq obtained from the above Vc{rY^'^^ are given in 
Fig. iroT Right). At low energy, the phase shift increases due to the attraction at medium distance, while 
at high energy it decreases as a consequence of the repulsive core at short distance. The shape of the 
scattering phase shift as a function of energy is qualitatively similar to but is much smaller than in 
magnitude the experimental one, plotted by the black solid line in Fig. [TOTRight). As discussed before, 
the pion mass dependence is so large for the scattering phase shift that full QCD simulations at physical 
pion mass are needed to reproduce the experimental behavior. 

It is noted here that ground state saturation has to be achieved to an accuracy of around 1 MeV, 
which is about 0.05 % of the total mass of the two-nucleon system, in order to determine the shift 
of the potential {E = k'^/rriN) at the same accuracy. The value of E has a strong influence on the 
value of the scattering length. Such a high precision is not yet attained in full QCD calculations, since 
significantly larger t and, accordingly, larger statistics are required. An alternative method to overcome 
this difficulty will be discussed in Sec. [7] . 
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Figure 11: (Left) 2+1 flavor QCD results for the spin-triplet central potential Vc{rY^'^^ from the orbital 
Af — T^ coupled channel at three values of the pion mass. (Right) The tensor potential Vrir) at three 
values of the pion mass. Taken from Ref. [35]. 



5 Hyperon Interactions 

Hyperon(y) potentials (hyperon-nucleon and hyperon-hyperon) serve as the starting point in studying 
the hyper-nuclei physics. Properties of these potentials can also determine structures in the core of 
neutron stars. In spite of their importance, only a limited knowledge for hyperon potentials is available 
so far, since experimental data such as scattering phase shifts are difficult to obtain, due to the short 
lifetime of hyperons. Therefore it is important to calculate hyperon potentials in lattice QCD by using 
the potential method. 

5.1 Quenched result for NE potentials 

Since all octet-baryons and decuplet Q are stable in the strong interaction, there are many hyperon 
potentials in 2 + 1 flavor QCD. The method for the NN potentials can be straightforwardly applied 
to the 1 = 1 NE channel, since pS is simply obtained from pn by replacing a d-quark in the neutron 
with the s-quark and it does not have strong decay into other channels. Unstable channels such as the 
7 = N'E, which can decay into AA via the strong interaction, will be discussed later. In addition, 
experimentally, not much information has been available on the iVS interaction except for a few studies: 
a recent report gives the upper limit of elastic and inelastic cross sections [36] while earlier publications 
suggest weak attractions of S— nuclear interactions [371 EHl ES]- The S— nucleus interactions will be 
soon studied as one of the day-one experiments at J- PARC [10] via {K~,K~^) reaction with a nuclear 
target. 

Ref. [H] gives the flrst quenched result for 1 = 1 N'E potentials. Lattice parameters are the same 
as for the quenched A^A^ potential. In addition to two values of the light quark mass, the quenched 
strange quark is introduced and is flxed to one value. The potential is calculated at {mT,,miy,m^) = 
(511(l)MeV,1300(4)MeV,1419(4)MeV) and (368(l)MeV, 1167(7)MeV, 1383(6)MeV), using the NBS 
wave function with the interpolation operators deflned by 



Paix) = eabciu''ix)C-f5d\x))u''^{x), S^(x) = eabciu''{x)C-f5S {x))s''^{x 



(47) 



Since both p and S° have (J, J^) = 1/2, 1/2, the pE^ system has 7=1 with the strangeness S = —2. 

The left (right) of Fig. [12] gives the (effective) central potential of the pS° system obtained from the 
L = Ai representation for the spin-singlet (triplet) at m^ = 511 MeV and 368 MeV. Potentials in the 
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Figure 12: (Left) The spin-singlet central potential for pS° obtained from the orbital Af channel at 
m^ ~ 368 MeV (circle) and m^ ~ 511 MeV (box). The central part of the OPEP {F/{F + D) = 0.36) in 
Eq. ( 148|) is also given by solid line. (Right) The spin-triplet effective central potential from the orbital 
Af channel at 771^^ ~ 368 MeV (triangle) and r7i,r — 511 MeV (diamond), together with the OPEP (solid 
line). Taken from Ref. |41j . 



1=1 NS system for both channels show a repulsive core at r < 0.5 fm surrounded by an attractive 
well, similar to the NN systems. In contrast to the A^A^ case, however, the repulsive core of the pS° 
potential in the spin-singlet channel is substantially stronger than in the triplet channel. The attraction 
in the medium to long distance region( 0.6 f m < r < 1.2 fm ) is similar in both channels. The height 
of the repulsive core increases as the light quark mass decreases, while a significant difference is not 
seen for the attraction in the medium to long distance within statistical errors. Potentials in Fig. [T21 
are weakly attractive on the whole in both spin channels at both pion masses, in spite of the repulsive 
core at short distance, though the attraction in the triplet is a little stronger than that in the singlet. 
The solid lines in Fig. [T2]are the one-pion exchange potential (OPEP), given by 



TATT 



;i - 2a) 



Att 3 



m. 



2m 



N 



(4J 



with 



m^,mN 



(368MeV, 1167MeV), where the pseudo-vector ttSS coupling /tthh is related to the 



/. 



NN- 



ttNN coupling as /ttsh = —fnNNi^ — 2a) with the parameter a = F/{F + D), and Qj^nn — jtti\i\ 2m ■ 
The empirical vales, a ~ 0.36 and g-^NN ji/^T^) — 14.0, are used for the plot. Unlike the A^A^ potential, the 
OPEP in the present case has opposite sign between the spin-singlet channel and spin-triplet channel. 
In addition, the absolute magnitude is smaller due to the factor 1 — 2a. No clear signature of the OPEP 
at long distance (r > 1.2 fm) is observed in Fig. [T2] within statistical errors. Furthermore, there is 
clear departure from the OPEP at medium distance (0.6 f m < r < 1.2 fm) in both channels. These 
observations may suggest an existence of state-independent attraction. 



5.2 Full and quenched QCD results for A^A potentials 

Spectroscopic studies of the A and S hypernuclei, carried out both experimentally and theoretically, 
suggest that the A-nucleus interaction is attractive while the S-nucleus interaction is repulsive. If this 
is the case, the A particle would be the first strange baryon instead of S^ to appear in the core of 
the neutron stars[l2]. It is therefore interesting and important to investigate the nature of the A^A 
interaction in lattice QCD, by calculating the A^A potential using the method of this report. Since the 
A is the lightest hyperon, the A^A potential can be calculated as in the case of A^A^ potentials. 
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Figure 13: (Left) The spin-singlet central potential for NA obtained from the orbital Af channel in 2+1 
flavor QCD at ?7i^ ~ 414 MeV (red) and 699 MeV (green). (Right) The spin-triplet central potential 
and the tensor potential for A^A obtained from the orbital Af — T2 coupled channel in 2+1 flavor QCD 
at TTi^ ^ 414 MeV (red and blue) and 699 MeV (green and magenta). Taken from Ref. |13] . 
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Figure 14: (Left) The spin-singlet central potential for A^A obtained from the orbital Af channel in 
quenched QCD at 171^^ ~ 407 MeV (red) and 512 MeV (green). (Right) The spin-triplet central potential 
and the tensor potential for A^A obtained from the orbital Af — T2 coupled channel in quenched QCD 
at m^ ~ 407 MeV (red and blue) and 512 MeV (green and magenta). Taken from Ref. |13]. 



In Ref. [13], the A^A potentials are calculated in both full and quenched QCD. The 2+1 flavor full 
QCD gauge conflgurations generated by the PACS-CS collaboration are employed for the calculations of 
the potentials on a 32'^ x 64 lattice at a = 0.091(1) fm, while in the quenched calculation, the potentials 
are obtained on a 32'^ x 48 lattice at a = 0142(1) fm. Numerical values for some hadron masses for 
these calculations are given in Table HI together with some lattice parameters. 

Fig. [T3l shows the A^A potentials obtained from 2+1 flavor QCD calculations as a function of r at 
mTT ^ 699 MeV and 414 MeV. The spin-singlet central potential obtained from the J = Ai channel is 
plotted on the left, while the spin-triplet central potential and the tensor potential obtained from the 
J = Ti channel are given on the right. The central potential multiplied by a volume factor {r^Vc{r)) 
is also shown in the left panel in addition to the Vc{t) itself in the right panel, in order to compare the 
strength of the repulsion between two quark masses. 

As can be seen in Fig. [131 the attractive well of the central potentials moves to the outer region as 
the light quark mass decreases, while the depth of these attractive pockets does not change so much. 
The present results show that the tensor force is weaker than the A^A^ case in Figj9l Moreover the quark 
mass dependence of the tensor force seems small. Both repulsive and attractive parts of the central 
potentials increase in magnitude as the hght quark mass decreases. 

Fig. [TH shows the A^A potentials in quenched QCD calculations at 771^^ ~ 512 MeV and 407 MeV. 
The central potential in the spin-singlet channel from J = Ai is on the left, while the central and 
the tensor potential in the triplet channel from J = Ti are on the right. Qualitative features of these 
potentials are more or less similar to those in full QCD: the attractive pocket of the central potentials 
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m^ 



nir 



niK 



niK* 



niN 



niA 



ms 



m-= 



2+1 flavor QCD on a 32^ x 64 lattice at a = 0.091(1) fm 
699.4(4) 1108(3) 786.8(4) 1159(2) 1572(6) 1632(4) 
567.9(6) 1000(4) 723.7(7) 1081(3) 1396(6) 1491(4) 
413.6(6) 902(3) 636.6(4) 1026(3) 1221(7) 1349(4) 
301(3) 845(10) 592(1) 980(6) 1079(12) 1248(15) 



1650(5) 1701(4) 

1519(5) 1599(4) 

1406(8) 1505(4) 

1308(13) 1432(7) 



quenched QCD on a 32^ x 48 lattice at a = 0.142(1) fm 
511.8(5) 862(3) 605.8(5) 898(1) 1297(6) 1344(6) 
463.6(6) 842(4) 586.3(5) 895(2) 1250(9) 1314(9) 
407(1) 820(3) 564.9(5) 886(3) 1205(13) 1269(9) 



1375(5) 1416(3) 
1351(6) 1404(4) 
1326(9) 1383(5) 



135 



770 



494 



892 



940 



1116 



1190 



1320 



Table 4: Hadron masses in units of MeV in lattice QCD simulations. The last line shows experimental 
values. Taken from Ref. 



moves to the longer distance region as the quark mass decreases while the quark mass dependence of 
the tensor potential is small. 

5.3 Flavor SU(3) limit 

In order to unravel the nature of the various channels in the hyperon interactions, it is more convenient 
to consider an idealized flavor SU(3) symmetric world, where u, d and s quarks are all degenerate with 
a common flnite mass. In the flavor SU(3) limit, one may capture essential features of the interaction, 
in particular, the short range force without contamination from the quark mass difference. Calculations 
in the SU(3) limit allow us to extract potentials for irreducible flavor multiplets: Potentials between 
asymptotic baryon states are obtained by the recombination of the multiplets with suitable Clebsh- 
Gordan coefficients. 

In the flavor SU(3) limit, the ground state baryon belongs to the flavor-octet with spin 1/2, and two 
baryon states with a given angular momentum are labeled by the irreducible representation of SU(3) as 



8®8 = 27©8©1©10©10©8, 



(49) 



symmetric 



anti— symmetric 



where "symmetric" and "anti-symmetric" stand for the symmetry under the exchange of the flavor 
for two baryons. For the system with orbital S-wave, the Pauli principle imposes 27, 8 and 1 to be 
spin-singlet {}Sq), while 10, 10 and 8 to be spin-triplet (^^i). 
The NBS wave function is deflned by 



^^'W(r) 



(0|(Si?)W(r)|W^,i? = 2,X), 



from which the corresponding (effective) central potential is given as 



rW(r 



1 

2^ 



V 



W,{X)i 



+ k' 



(50) 



(51) 



with k"^ = W"^ /A — m\ and yU = mBJ'^-, where m^ is the common octet baryon mass and {BB)'^-^^(j) = 
Y.ij CijBi{yi + r, 0)5j(x, 0) with X = 27, 8^, 1,10, 10, 8a. Two baryon operators BB^^^ in the flavor 
basis are given in terms of the particle basis in Appendix |Bl Potentials among octet baryons, both 
the diagonal part (-B1-B2 — > B1B2) and the off-diagonal part (-B1-B2 — ^ B3B4) are obtained by suitable 
combination of V^^\r). 
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Figure 15: The BB potentials in 27 (Left) and lO(Right) representations from the orbital Af channel 
in the flavor SU(3) limit, extracted from the lattice QCD simulation at m.^ = 1014 MeV (red bars) and 
m.^ = 835 MeV (green crosses). Taken from Ref. [15] . 



flavor multiplet 


baryon pair (isospin) 


27 

8. 
1 


{NN}(I=1), {NS}(I=3/2), {SS}(I=2), 

{SS}(I=3/2), {SS}(I=1) 

none 

none 


10* 
10 

8a 


NN (1=0), [SH (1=3/2) 
[NS] (1=3/2), [HS](I=0) 
[NH](I=0) 



Table 5: Baryon pairs in an irreducible flavor SU(3) representation, where {BB'} and [BB'] denotes 
BB' + B'B and BB' - B'B, respectively. 



Using the 3 flavor full QCD gauge configuration generated by CP-PACS/JLQCD Collaborations 
on a 16^ X 32 lattice at a ~ 0.12 fm where light and strange quark masses have same values[Hj, 
the (effective) central potentials are calculated |15] at {mps.rnB) = (1014(l)MeV, 2026(3)MeV) and 
(845(l)Mev, 1752(3)MeV), where mps and m^ denote the octet pseudo-scalar (PS) meson mass and 
the octet baryon mass, respectively. 

Figs. [I5] and [I6] give the six baryon-baryon {BB) potentials in the flavor basis, where red (green) 
data correspond to the mps = 1014 MeV (835 MeV). The left panels show central potentials for the 
spin-singlet channel from the J = Ai state, while the right panels give effective central potentials for 
the spin-triplet channel from the J = Ti state. Note that some of octet-baryon pairs exclusively belong 
to an irreducible representation of SU(3) as shown in Table O For example, symmetric A^A^ belongs to 
27, which therefore can be considered as the NN spin-singlet potential in the flavor SU(3) symmetric 
limit. Similarly V'^^^\ V^'^^'^ and V^(*") can be considered as some BB potentials of the particle basis in 
the SU(3) symmetric limit, while V^^"* and y(*=) are always mixtures of different BB potentials in the 
particle basis. 

Fig. [T5] shows V^'^'^\r) and V'^^^^r), which correspond to spin-singlet and spin-triplet NN potentials, 
respectively. Both have a repulsive core at short distance with an attractive pocket around 0.6 fm. These 
qualtative features are consistent with the previous results found for the NN potential in both quenched 
and full QCD. The upper-right panel of Fig. [16] shows that V^^^\r) has a stronger repulsive core and 
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Figure 16: The BB potentials in 8^ (Upper-Left), 10(Upper- Right), l(Lower-Left) and 8^ (Lower- Right) 
from the orbital Af channel in the flavor SU(3) limit, extracted from the lattice QCD simulation at 
m^r = 1014 MeV (red bars) and m^ = 835 MeV (green crosses). Taken from Ref. [iS] . 



a weaker attractive pocket than V^^'^'^^\r). Furthermore \/(*'')(r) in the upper-left panel of Fig. [TBI has 
a very strong repulsive core among all 6 channels, while \/(^")(r) in the lower-right panel has a very 
weak repulsive core. In contrast to all other cases, V^'^\r) shows attraction instead of repulsion at all 
distances, as shown in the lower-left panel. 

Above features are consistent with what has been observed in a phenomenological quark model|l6]. 
In particular, the potential in the 8^ channel in this quark model becomes strongly repulsive at short 
distance since the six quarks cannot occupy the same orbital state due to the Pauli exclusion for quarks. 
On the other hand, the potential in the 1 channel does not suffer from the quark Pauli exclusion and 
can become attractive due to the short-range gluon exchange. Such agreements between the lattice 
data and the phenomenological model suggest that the quark Pauli exclusion plays an essential role for 
the repulsive core in BB systems. 

The BB potentials in the baryon basis can be obtained from those in the SU(3) basis by a unitary 
rotation as 

V^Ar) = T.U.xV^''\r)ul,^ (52) 

X 

where [/ is a unitary matrix which rotates the flavor basis |X) to the baryon basis \i), i.e. \i) = Uix\X). 
The explicit forms of the unitary matrix U in terms of the CG coefficients are found in Appendix [Bl 

In Fig. [T71 as characteristic examples, let us show the spin-singlet potentials for S=— 2, 1=0 channel 
determined from the orbital Af representation at m„ = 835 MeV. To obtain V^j(r), the potentials in 
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Figure 17: BB potentials in baryon basis for the S=— 2, 1=0, ^5*0 sector. Three diagonal (off-diagonal) 
potentials are shown in left (right) panel. Taken from Ref. [45j . 



the SU(3) basis are fitted by the following form with five parameters 61,2,3,4,5 



V{r) = bie-''''' +h{l-e-'^'' ) 



64 r 



2, / e' 



-fesr ' 



(53) 



Then the right hand side of Eq. ( |52l) is used to obtain the potentials in the baryon basis. The left panel 
of Fig. [T7] shows the diagonal part of the potentials. The strong repulsion in the 8^ channel is reflected 
most in the SS(I=0) potential due to its largest CG coefficient among three channels. The strong 
attraction in the 1 channel is reflected most in the iVS(I=0) potential due to its largest CG coefficient. 
Nevertheless, all three diagonal potentials have a repulsive core originating from the 8^ component. The 
right panel of Fig. [17] shows the off-diagonal parts of the potentials which are comparable in magnitude 
to the diagonal ones. Since the off-diagonal parts are not negligible in the baryon basis, fully coupled 
channel analysis is necessary to study observables. A similar situation holds even in (2-M)-fiavors where 
the strange quark is heavier than up and down quarks: The SU(3) basis with approximately diagonal 
potentials is useful for obtaining essential features of the BB interactions, while the baryon basis with 
substantial magnitude of the off-diagonal potentials is necessary for practical applications. 

Other potentials in the baryon basis are given in Ref. [15]. Since the 8^ state does not couple to the 
spin-triplet channel, the repulsive cores in the spin-triplet channel are relatively small. The off-diagonal 
potentials are not generally small: For example, the NA-NT, potential in the spin-triplet channel 
is comparable in magnitude at short distances with the diagonal NA-NA and NH-NTj potentials. 
Although all quark masses of 3 flavors are degenerate and rather heavy in these simulations, the coupled 
channel potentials in the baryon basis may give useful hints for the behavior of hyperons (A, S and S) 
in hyper- nuclei and in neutron stars |47l |48] . 

The flavor singlet channel has attraction for all distances, which might produce the bound state, the 
if-dibaryon, in this channel. The present data, however, are not sufficient to make a definite conclusion 
on the if-dibaryon, since a single lattice with small extension L ~ 2 fm is employed. In order to 
investigate whether the if-dibaryon exists or not in the flavor SU(3) limit, data on several different 
volumes are needed. Such a study on the if-dibaryon will be discussed in Sec. [71 

In order to extend the study in the flavor SU(3) limit to the real world where the strange quark is 
much heavier than light quarks, the potential method used so far has to be extended to more general 
cases, which will also be considered in Sec. [71 
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6 Origin of repulsive core 

As seen in the previous sections, lattice QCD calculations show that the A^A^ potential defined through 
the NBS wave function has not only the attraction at medium to long distance that has long been well 
understood in terms of pion and other heavier meson exchanges, but also a characteristic repulsive core 
at short distance, whose origin is still theoretically unclear. Furthermore, the BB potentials in the fiavor 
SU(3) limit show several different behaviors at short distance: some have a stronger /weaker repulsive 
core than A^A^ while the singlet has an attractive core. In this section, recent attempts [IHl EDI E] to 
theoretically understand the short distance behavior of the potential in terms of the operator product 
expansion (OPE) is explained. 

6.1 OPE and repulsive core 

Let us first explain the basic idea. We consider the equal time NBS wave function defined by 

^ABir) = {0\T{OAir/2,0)OBi-r/2,0)}\E), (54) 

where \E) is some eigenstate of a certain system with total energy E, and Oa, Ob are some operators 
of this system. (We suppress other quantum numbers of the state \E) for simplicity.) The OPE reads 

0A{r/2, 0)OB{-r/2, 0)c^J2 D%{r)Oc{0, 0), (55) 

c 

Suppose that the coefficient function of the OPE behaves in the small r(= \r\) limit as 

D^sir) ^ r^^{-\ogrf^fci9,(P), (56) 

where 6, are the angles of r, the NBS wave function becomes 

^IbM ^ Y.'^^^i-XogrY^ fc{e,<P)Dc{E), (57) 

c 

where 

Dc{E) = (0|Oc(0,0)|E). (58) 

The potential at short distances can be calculated from this expression. For example, in the case of the 
Ising field theory in two dimensions, the OPE for the spin field a is given by 

a(a;,0)or(0,0) ~G'(r)l + cr3/^Oi(0) + ---, r = |x|, (59) 

where Oi{x) (=: ipiplx) : in terms of free fermion fields) is an operator of dimension one. This leads to 

(^(r,E) ~r=^/^D(E) + 0(r^/^), D{E) = c{0\Oi{0)\E), (60) 



where \E) is a two-particle state with energy E = 2\^W+m?. From this expression the potential 
becomes 

mip{r, E) 16 mr^ 

in the r — )■ limit. The OPE predicts not only the r~^ behavior of the potential at short distance but 
also its coefficient —3/16. Furthermore the potential at short distance does not depend on the energy 
of the state in this example [30] 152] . 

In QCD the dominant terms at short distance have ac = 0. Among these terms, we assume that 
C has the largest contribution such that (3c > Pc' for VC" ^ C . Since such dominant operators with 
ac = mainly couple to the zero angular momentum (L = 0) state, let us consider the NBS wave 
function with L = 0. Applying V^ to this wave function, we obtain the following classification of the 
short distance behavior of the potential. 
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(1) /^c* 7^ 0: The potential at short distance is energy independent and becomes 

\/(r)^ /^ -, (62) 

which is attractive for Pc > Q and repulsive for /3c < 0. 

(2) Pc = 0: In this case the potential becomes 

where Pc' < is the second largest exponent. The sign of the potential at short distance depends 
on the sign of Dc>{E)/Dc{E). 

On the lattice, we do not expect divergence at r = due to lattice artifacts at short distance. The 
above classification holds at a ^ r ^ I/Aqcd, while the potential becomes finite even at r = on the 
lattice. 

Since QCD is an asymptotic free theory, the 1-loop calculation for anomalous dimensions becomes 
exact at short distance. The OPE in QCD is written as 

OA{y/2)OB{-yl2) = Y.D%{r,g,m,fx)Oc{0) (64) 

c 

where g (rra) is the renormalized coupling constant (quark mass) at scale /i. In the limit that r = \y\ = 
e~*-R — > ( t — > oo with fixed R), the renormalization group analysis leads to 



\imD^^{r,g,m,fi) = (-2/3«/logr)^^f V(2^«)/)C^(^^ q, 0,/i), (65) 



where (3^^^ = ^r (11 ] is the QCD beta- function at 1-loop, and 

167r^ V 3 / 

c.(i) (1) (1) (1) — ^ /««\ 

lAB =rc -Ta -Tb =-rr^l- (66) 

Here 'jx is the 1-loop anomalous dimension of the operator Ox- An appearance of Z)^g(i?, 0, 0,/i) 

on the right-hand side tells us that it is enough to know the OPE only at tree level. From the above 

c,(i) 
expression, f3c is given by Pc = l^J)- 

6.2 Two flavor case 

We first consider the OPE for Nf = 2 QCD|i49j. The 1-loop calculations show that the largest value of 
/3c is always zero for both spin-singlet and spin-triplet channels for Nf = 2 QCD and that the second 
largest value of /3c' is given by 

f^'^' = - 33 - 2N/ ^'^'^~33-2N/ ^^^^ 

where S* = 0, 1 denotes the total spin. This corresponds to the case (2) in the previous subsection. 
Therefore the OPE and renormalization group analysis in QCD predicts the universal functional form 
of the A^A^ central potential at short distance as 



Dc'{E)-(3i„{-\ogr[ 



P^.,-1 



Dc{E) ruNT^ 
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X 


27 


8. 


1 


10 


10 


8a 


Non-relativistic op. 




yes 


6 
no 


12 

yes 




yes 




yes 


4 
yes 



Table 6: The largest value of /3c in unite of 1/(33 — 2Nf) of 3-flavor QCD for each representation. 
The last line indicates that the operator corresponding to the largest value of Pc exists or not in the 
non-relativistic limit. 



which is a little weaker than a 1/r^ singularity, while for the tensor potential we have 

Vrir) ~ (69) 

at the 1-loop order. 

The OPE, however, can not tell whether the potential at short distance is repulsive or attractive, 
which is determined by the sign of the coefficient. If Dx{E) and Dy{E) are evaluated by the non- 
relativistic quark model at the leading order, we obtain 

For both cases, the ratio has positive sign, which gives repulsion at short distance, the repulsive core. 

6.3 Extension to three flavors 

The above calculation has been extended to Nf = 3 QCD [51]. In the 3-flavor case, some channel may 
become attractive at short distance since the Pauli exclusion principle is less significant than in the 2- 
fiavor case. Indeed the lattice QCD calculations in the fiavor SU(3) limit shows the attractive potential 
for the singlet channel, as seen in the previous section. 

The largest value of /3c is given for each X representation of the fiavor SU(3) in unit of 1/(33 — 2Nf) 
in the table [6l where we define 

"- = 33^- P^) 

From the table, we observe that the largest value of /3c is zero in the 27, 10 and 10 channels. This 
is consistent with the nucleon case in the previous subsection, which belong to 27 (spin-singlet) and 
lO(spin-triplet). These three channels correspond to the case (2), so that the potentials are given by 
eq. ( 163|) . On the other hand, the largest value of /3c becomes positive in the 8^, 8^ and 1 channels, 
which correspond to the case (1). Therefore the (effective) central potential becomes attractive at short 
distance as 

(X) 7^^^ 1 
V^ '(r) ~ - - ('72'! 

^ ^ ' {33-2Nf)mBr^{-logry ^ ' 

where m^ is the octet baryon mass. 

The attractive core of the potential in the flavor singlet channel agrees with the behavior of the 
potential found for the numerical simulation of lattice QCD in the previous section, while for other two 
channels, 8s and 8a, the prediction by the OPE disagrees with the lattice QCD results: The potential in 
the 8s channel is most repulsive among 6 channels and the potential in the 8a channel still has a repulsive 
core, which is however weaker than others. The disagreement between the OPE and the lattice QCD 
result for the 8s channel may be understood by the fact that no local 6 quark operator exists for this 
channel in the non-relativistic limit, as shown in the table [6j the 8s operator with the largest positive 
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Pc has a very small coefficient at low energy, so that the other operators with zero or negative Pc ma-y 
still dominate at distance scales comparable to the lattice spacing a = 0.1 — 0.2 fm. For the 8^ case, 
the weakest repulsive core in the lattice QCD simulation suggests that the attraction from the leading 
operator with the positive /3c may be cancelled by other contributions from the sub-leading operators 
with zero or negative f3c at the distance scale comparable to the lattice spacing of the simulations. It 
is therefore important to confirm the prediction from the OPE, by investigating the behavior of the 
repulsive core for each channel in the fiavor SU(3) limit at finer lattice spacings and hopefully in the 
continuum limit. 

7 Extensions 

In this section, recent extensions of the potential method are considered. 

7.1 Inelastic scattering 

The potential method discussed so far is shown to be quite successful in order to describe elastic hadron 
interactions. Hadron interactions in general, however, lead to inelastic scatterings as the total energy of 
the system increases. In order to extract hadron interactions which describe such inelastic scatterings 
from lattice QCD, an extension of the potential method is considered in this subsection. 

Let us first discuss the case oi A + B — t- C + D scattering where A,B,C,D represent some 1-particle 
states. This is a simplified version of the octet baryon scattering in the strangeness S = —2 and isospin 
/ = channel, where AA, NE and SS appear as asymptotic states of the strong interaction if the total 
energy is larger than 2ms. 

We here assume m^ + m^ < nic + m£, < W, where W = E^ + E^^ is the total energy of the system. 



and E^ = \lm\ + A; . In this situation, the QCD eigen-state with the quantum numbers of the AB 
state and center of mass energy W is expressed in general as 

\W) = cab\AB,W) + ccd\CD,W) + --- (73) 

\AB,W) = |A,fc)i,®|fi,-fc)i„ \CD,W) = \C,q),^®\D,-q),^, (74) 

where W = E^ + E^ = E^ + Ej^. We define the following NBS wave functions, 

ipAB{r,k)e-'^' = {0\T{ipA{x + r,t)ipB{x,t)}\W), (75) 

cpcDir,q)e-'^' = {0\T{ipcix + r,t)ipnix,t)}\W). (76) 

Using the partial wave decomposition such thacl 

ipxYir,k) = ATTj2i^^xYir,k)Yi„,inr)Yi^ink), (77) 

l,m 

the NBS wave function of the 2-channel system behaves for large r as 

f^AB{r,k)\ _ ( Ji{kr) \ ( cab \ ^ ( ni{kr) + iji{kr) A 

V fcDir^q) J \0 Jiiqr) J \ Ccd J \0 ni{qr) + iji{qr) J 

"" ^^^n e^^?(^)sin5f(iy) J^ ^^\ccn)' ^^^^ 

^ ' \ sva.6{W) cos6'(iy) / ^ ' 



*Here we ignore spins for simplicity. 
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where Sl{W) is the scattering phase shift , whereas 6{W) is the mixing angle. This expression shows 
that the NBS wave functions for large r agree with scattering waves described by two scattering phases 
Sl{W) {i = 1, 2) and one mixing angle 9(W). Because of this property, these wave functions satisfy 



(V^ + e)^AB{r, k) = 0, (V' + q^)^cD{r, q) 







^0) 



for r — )■ oo. 

Let us now consider QCD in the finite volume V. In the finite volume, \AB, W) and \CD, W) are 
no longer eigen-states of the hamiltonian. True eigenvalues are shifted from W to Wi = W + 0{V~^) 
{i = 1,2). By diagonalization method in lattice QCD simulations, it is relatively easy to determine 
Wi and W2- With these values Liischer's finite volume formula gives two conditions, which, however, 
are insufficient to determine three observables, 6}, 6f and 9. (See [53l EU [55] for recent proposals to 
overcome this difficulty.) An alternative approach to extract three observables. Si, Sf and 6, has been 
proposed in lattice QCD through the above NBS wave functions[56l [57]. We consider the NBS wave 
functions at two different values of energy, Wi and W2, in the finite volume: 



(^CD(r,g,)e-^»* 



{0\T{^A{x + r,t)^B{x,t)}\Wi) 
{0\T{ipcix + r,t)ipDix,t)}\Wi), % 



1,2. 



(81) 
f82) 



We then define the coupled channel non-local potentials from the coupled channel Schrodinger equation 
as 



2/i^B 



E, 



ipAB{x, ki) = J d^y Uab,ab{x; y) (fAB{y, ki) + / d^y Uab,cd{x; y) (ficniy, Qt 



^HCD 



-H^ 



ipCD{x,ki) 



d^y UcD,AB{x]y) ^AB{y,ki) + / d^y UcD,CD{x]y) ^cD{y,q^ 



for z = 1, 2. As before the velocity expansion is introduced as 

UxY,vz{x; y) = Vxy,vz{x, V)5^ix - y) = [Vxy,vz{x) + 0(V)] 5\x - y) 
and at the leading order of the expansion, we have 



(83) 

I 

(84) 
(85) 



KAB{x,ki) = 



KcD{x,qi) 



'^jJ.AB 



Hn 



(Pab{x, ki) = Vab,ab{x) (fABix, ki) + Vab,cd{x) (Pcd{x, gj 



2/icz) 



-Ho 



ipcD{x,ki) 



Vcd,ab{x) (Pab{x, ki) + Vcd,cd{x) ipcoix, qj 



^6) 



^7) 



These equations for i = 1, 2 can be solved as 

f Vab,ab{x) Vab,cd{x) 
V Vcd,ab{x) Vcd,cd{x) 



X 



KABix,ki) KAB{x,k2) 

KcD{x,qi) KcD{x,q2) 

ipAB{x,ki) ipAB{x,k2) 

ifCDix^q^) ipcD^x^q^) 



Once we obtain the coupled channel local potentials Vxy,vz{x), we solve the coupled channel 
Schrodinger equation in infinite volume with some appropriate boundary condition such that the in- 
coming wave has a definite I and consists of the AB state only , in order to extract three observables 
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for each (i {5}{W), ^fiW) s-nd 0{W)) at all values of W . Of course, since Vxy,vz is the leading order 
approximation in the velocity expansion of UxYyz{.x-^y)i results for three observables 5}{W), 5f(W) 
and 9{W) aXW ^ Wi, W2 are also approximate ones and might be different from the exact values. By 
performing an additional extraction of Vxy,vz{x) at (1^3,^4) 7^ (Wi,W2), we can test how good the 
leading order approximation is. 

The method considered above can be generalized to inelastic scattering where a number of particles 
is not conserved. For illustration, let us consider the scattering A + B — > A + B and A + B — )■ A + B + C 
where the total energy W satisfies rriA + thb + rnc <W < niA + trb + 2mc. 

The following NBS wave functions at the center of mass system are used: 



where 






(0|T{v9A(r4 



x,t)^B{r,t)}\W) 

^ ^ + ^if^, t\ ^^^r + y, t)Mr, t)}\W), 
mc J 



\W) = Ci |fc)in8) I -fc)in + C2 |q^)i 






Qxf^BcX 



VflB 



(89) 
(90) 



(91) 



with 



W 



k' 



m, 



k^ + m| 



ql + m\ + Mq - 



mc 



l^^ + 'ki^y + ml 

^ m,B 



(92) 



and l/fiAB = i/m^B + 1/m.c- Here 7/ = r^ — re is a relative coordinate between B and C with the 
reduced mass fiBC, while x = va — Rbc is the one between A and the center of mass of B and C with 
Rbc = (mBTB + mcrc)/{mB + mc). 

We define the non-local potential from the coupled channel equations as 



K 



w 

AB 



X] 



KABci^^y) 



where 



k' 



'^fJ'AB 



H, 



AB 



V^Ab(*) = j ^^^ UaB,Ab{x] z) ifABi^) 



+ I £z d^w Uab,abc{x; z, w) (^ABci^^ ^) 



(93) 



+ 



-H, 



A,BC 



H, 



BC 



X 



'^fJ'A,BC '^f^BC 

VabI^) + J d^ z d^w Uabc,abc{x, y; z, w) if^Bci^^ '^ 



^ABci^^ y) = J d^z Uabc,ab{x, y; z) 



H, 



AB 



V2 

^ X TT 

7, > -"0 

^IJ'AB 



A,BC 



'^f^A,BC 



H, 



BC 



^ y 



'2fJ,BC 



with another reduced mass defined by l//iA,BC = ^/m.A + 1/(^b + m^c)- 
We consider the following velocity expansions 



UAB,ABix;z) 

Uab,abc{x;z,w) 

UABc,AB{x,y;z) 

UABC,ABcix,y;z,w) 



[VabMx) + 0{W)]S%x-z) 
[Vab,abc{x, w) + 0(V")] 6^{x - z) 
[VAB,ABc{x,y) + 0{V'')]5\x-z) 
[Vabc,abc{x, y) + 0(V^ V^)] 5\x - z)5\y - w), 



(94) 



(95) 
(96) 
(97) 
(98) 
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Set 1 
Set 2 
Set 3 



-'■ " conj 

700 
800 
800 



m. 



niK 



lllN 



niA 



m^ 



m-^ 



875(1) 916(1) 1806(3) 1835(3) 1841(3) 1867(2) 
749(1) 828(1) 1616(3) 1671(2) 1685(2) 1734(2) 
661(1) 768(1) 1482(3) 1557(3) 1576(3) 1640(3) 



Table 7: Hadron masses in units of [MeV] and number of configurations for each set. 



where the hermiticity of the non-local potentials gives Vab,abc{x, y) = Vabc,ab{x, v)- 
At the leading order of the velocity expansions, the coupled channel equations become 



K 



w 

AS 



X] 



Vabab{.x)'{>^b{x) + I d^w Vab,abc{x,w)ip^ 



^^(jyX,W) 



^W 



w 



^ABci^^ y) = yABc,AB{x, y)^^^(a;) + Vabc,abc{x, y)ipABci^, y). 



(99) 
(100) 



By considering two values of energy such that W = Wi, W2, we can determine Vabc,ab and Vabc,abc 
from the second equation as 



( VABC,ABix,y) VABC,ABcix,y) 



K^^ci^,y) KTBc{^,y 



W2 



X 






x] 



^ 



nhci^.y) nici^.y 



W2, 

AB^ 
W2 



X 



(101) 



Using the hermiticity relation VAB,ABc{x,y) = VABC,ABix,y), we can extract Vab,ab from the first 
equation as 



VAB,ABix) 



^ 



W 
AB 



X 



KJb{^) - j d^w Vabc,ab{x,w)'^^jj(.{x,w) 



(102) 



for W = Wi, W2. A difference of Vab,ab{x) between two estimates at Wi and W2 gives an estimate for 
higher order contributions in the velocity expansions. 

Once we obtain Vab,ab, Vab,abc = Vabc,ab and Vabc,abc, we can solve the coupled channel 
Schrodinger equations in the infinite volume, in order to extract physical observables. As W increases 
and becomes larger than mA+iTLB+nmc, the inelastic scattering A+B — )■ A + B+nC becomes possible. 
As in the case oi A + B — )■ A + B + C in the above, we can define the coupled channel potentials including 
this channel, though calculations of the NBS wave functions for multi-hadron operators become more 
and more difficult in practice. 



7.2 Coupled channels with S = —2 and 7 = 

As an application of the method in the previous subsection, let us consider BB potentials for the S = —2 
and / = channel, which consist of the AA, NS and SS components in terms of low-lying octet baryons. 
Mass differences of these components are quite small such that 2mx = 2232 Mev, rriN + mj^ = 2257 
MeV and 2ms = 2386 MeV. Using diagonalized source operators, the NBS wave functions at three 
different values of energy. 



^w. 



Vabv^ ^ab)^ 



'Wd 



{Q\T{^A{r + x,t)^B{r,t)}\Wi) 



(103) 



for i = 0, 1, 2, are extracted in lattice QCD simulations, where AB = AA, NS and SS, and k^^ satisfies 



W^ = J{k\^y + m\ + J{k\^ 



m 
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Figure 18: The coupled channel potential matrix from the NBS wave function for Set 1. The vertical 
axis is the potential strength in units of [MeV], while the horizontal axis is the relative distance between 
two baryons in units of [fm]. Taken from Ref. 



Using the notation 



K'ABir) 



2/i 



AB 



(V^ + (fcAB)>TB(^,fcAB) 



where l/fiAs = l/'^A + l/'^s? the coupled channel 3x3 potential matrix is given by 



(105) 



(106) 



Here the last factor is the inverse of the 3x3 matrix (^^^^{r, k^^^) with indices i and CD. 

Gauge configurations generated on a 16^ x 32 lattice at a ~ 0.12 fm ( therefore L ~ 1.9 fm) in 2+1- 
flavor full QCD simulation are employed to calculate the coupled channel potentials at three different 
values of the light quark mass with the fixed bare strange quark mass [58] • Quark propagators are 
calculated with the spatial wall source at to with the Dirichlet boundary condition in time at t = to + 16. 
The wall source is placed at 16 different time slices on each gauge configuration, in order to enhance 
signals, together with the average over forward and backward propagations in time. Corresponding 
hadron masses and number of gauge configurations are given in table [3 

The coupled channel potential matrix Vab,cd from the NBS wave function for Set 1 is shown 
in Figure [13 The flavor dependence of the height of the repulsive core at short distance region is 
observed. In particular, the SS potential has the strongest repulsive core of these three channels. It is 
interesting to see that off-diagonal parts of the potential matrix roughly satisfy the hermiticity relation 
Vab,cd = VcD,AB within statistical errors. In addition the off-diagonal parts are similar in magnitude 
for Vaa,se and Vnb.,j]t. with the diagonal parts, but Vaa,ne. is much smaller. 
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Figure 19: Transition potentials in the flavor SU(3) IR basis. Red, blue and green symbols correspond 
to results of Setl, Set2 and Set3, respectively. The result of the flavor SU(3) symmetric limit at the 
same strange quark mass is also plotted with brown symbols [15]. Taken from Ref. [58] . 



In order to compare the results of the potential matrix calculated in three configuration sets, the 
potentials from the particle basis are transformed to those in the flavor SU(3) irreducible representation 
(IR) basis as 



V 



IR 



^1,1 ^1,8 ^1,27 

U^VU = I ^8,1 ^8,8 ^8,27 

^7,1 ^27,8 ^27,27 



(107) 



where t/ is a unitary transformation matrix whose explicit form is given in Appendix |Bl The potential 
matrix in the IR basis is convenient and a good measure of SU(3) breaking effects by comparing three 
configuration sets since it is diagonal in the SU(3) symmetric limit. 

In Figure [191 the results of the potential matrix in the IR basis are compared between different 
configuration sets, together with the one in the flavor SU(3) symmetric limit. As the pion mass decreases, 
the repulsive core in the ^27,27 potential increases. The 1^1^27 and ^§,27 transition potentials are consistent 
with zero within statistical errors. On the other hand, it is noteworthy that the flavor SU(3) symmetry 
breaking effect becomes manifest in the Vi^s transition potential. 



7.3 Time dependent method 

One of the practical difficulties to extract the NBS wave function and the potential from the correlation 
function Eg. (12^ is to achieve the ground state saturation in numerical simulations at large but finite 
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t — to with reasonably small statistical errors. While the stability of the potential against t — to has been 
confirmed within statistical errors in numerical simulations reviewed in the report, the determination of 
W for the ground state suffers from systematic errors due to contaminations of possible excited states, 
which can be seen as follows. There exist 3 different methods to determine W. The most well-known 
method is to determine W from the t — to dependence of the correlation function eg. ( 12^ summed over 
r to pick up the zero momentum state. On the other hand, one can determine k^ of W by fitting the 
r dependence of the NBS wave function with its expected asymptotic behavior at large r or by reading 
off the constant shift of the Laplacian part of the potential from zero at large r. Although the latter 
two methods usually give consistent results within statistical errors, the first method (the t dependence 
method) sometimes leads to a result different from those determined by the latter two (called the r 
dependent method together) at the value of t — to usually employed in numerical simulations. Although, 
in principle, the increase of t — to is needed in order to see an agreement between t and r dependence 
methods, it is difficult in practice due to larger statistical errors at larger t — to for the two-baryon 
system. 

In order to overcome this practical difficulty, the method to extract the potential from the NBS 
wave function has been modified as follows. Let us consider the correlation function Eq. fl2^ again: 



F{r,t)= Yl ^11/0^^(^)6-^* + 0(e-^*'^*). (108) 

W<Wt^ 

If t is large enough so that contributions from 0(e-^"^*) terms can be neglecteqj, we have 

HoF{r, t)^J2^w [ d^r'[EwS^^\r - r') - U{r, r')]ip^ {r')e-^' (109) 

where Ew = kyy/{2n) = {W^ — 4m^)/(4m7v) with /i = 171^/2. By using the non-relativistic approxi- 



mation that W = 2^JkiY + mj^ = 2mjv + k^y/m]\f + 0{k^r/m%), 

F{r, t) = - [ d^r'U{r, r')F{r, t) ~ -V^^{r)F{r, t) (110) 



rr d ^ 

Ho + — + 2mN 
dt 



where the velocity expansion is introduced in the last line and higher other than the leading order terms 
are then omitted. The leading order potential is therefore given by 



V^'^ir) 



Ho + i + 2mN F{r,t) 



F{r,t) 
or 



V^^(r) 



R{r,t) 



:iiii 



(112) 



R{r,t) 

where R{r,t) = F(r, t)/e-^™^*. Here it is assumed that 0(e-^"'*) contributions can be neglected at 
large t. The non-relativistic formula for V^^{r) above can be easily generalized to the case that masses 
of two particles are different, by the replacement that R{r,t) = F{r,t)/e~^"^^~^"^'^^^. Note also that 
the potential extracted in this method automatically satisfies that V^^*-*(r) — )■ as r — ;■ without the 
constant shift. This property may be used to check whether this extraction works correctly or not. 

On the lattice, the t derivative should be approximated by the t difference. In practice, one may 
adopt a particular method for the t difference, in order to reduce statistical as well as systematic errors 
for V^^{r). 



^This limitation for t can be loosened if the coupled channel potentials are introduced as in the previous subsections. 
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The non-relativistic approximation can be removed by using the second order derivative in t as 

F{r,t) 



V^^ir] 



^ 1 £ 
-Ho + — -niN 



4m AT dt"^ 



F{r,t) 



(113) 



as long as 0(e~^**) contributions are neghgibly small. For this method to apply, two particles should 
have the same mass. Statistical errors of the second order difference on the lattice must be kept small 
in numerical simulations. 

One may introduce a more general correlation function as 



F{x,y,t)= d'x^d'y,{0\T{N{x, + x,t)N{x^,t)}T{N{y, + y,0)N{y„0)}\0). 



Using this new quantity, we have 

H ^ ^' 

-no 



+ niN 



Arrijsf dt"^ 
from which the non-local potential is extracted as 



R{x, y,t) = -Jd^ zU{x, z)F{z, y, t). 



(114) 



;ii5) 



Uix,y) 



d-'z 



1 rf2 

-Hq + — -rriN 



F{x,z,t)-F-\z,y,t). 



4m iv dt^ 
Here F^^{x, y, t) is the approximated inverse of the hermitian operator F{x, y, t), defined by 

F~\x,y,t)= J2 Y7TT^n(a3,t)t;l(i/,t) 
x„^o ^"I'-J 



;ii6) 



(117) 



where A„(t) and v{x,t) are an eigenvalue and a corresponding eigenvector of F{x,y,t), respectively, 
and zero eigenvalues are removed in the summation. Since the modified potential 



U{x,y) = U{x,y) + ^ CnVn{x,t)vl{y,t) 

A„=0 



:ii8) 



also satisfies the same Schrodinger equation V{c„}, the non-local potential is NOT unique, and U{x, y) 
is scheme dependent, as discussed before. 

7.4 Bound H dibaryon in flavor SU(3) limit 

As an application of the method in the previous subsection, let us consider the singlet potential in the 
flavor SU(3) limit in order to investigate whether the bound H dibaryon exists or not in this case. 
At the leading order of the velocity expansion, the central potential is defined in this method by 



l^W(,) 






R{r,t-to) 



i?(r,t-to) 



(119) 



which is calculated on 16^ x 32, 24^ x 32 and 32^ x 32 lattices at a = 0.121(2) fm and three values 
of the quark mass, where the PS meson mass and the octet baryon mass are given by (mps,'"^^) = 
(1015(l)MeV,2030(2)MeV), (837(l)MeV, 1748(l)MeV) and (673(l)MeV, 1485(2)MeV) on a 32^ x 32 
lattice l59j. 
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MeV and (t - to)/a = 10. Taken from Ref. 
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Figure 21: The flavor- singlet potential V^^\r) at {t-to)/a = 10. (Left) Results for L = 1.94, 2.90, 3.87 
fm at mps = 1015 MeV. (Right) Resuhs for L = 3.87 fm at nip, = 1015, 837, 673 MeV. Taken from 
Ref. ESI. 



To check the qualitative consistency with previous results, the central potential in the 27-plet channel 
is plotted in Fig. |20] obtained in three different lattice volumes with L = 1.94, 2.90, 3.87 fm at rrips = 
1015 MeV and {t — to)/a = 10. This is the case corresponding to the spin-singlet NN potential. 
Compared with statistical errors, the L dependence is found to be negligible. The t dependence is also 
small as long as [t — tQ)/a > 9. As expected, the potential approaches zero automatically for large r. 
The figure shows a repulsive core at short distance surrounded by an attractive well at medium and long 
distances, which is qualitatively consistent with previous results in quenched and full QCD simulations. 

Shown in Fig. [2T](Left) and Fig. [2T](Right) are the volume and the quark mass dependences of the 
central potential in the flavor- singlet channel Vq (r), respectively, at {t—tQ)/a = 10 where the potentials 
do not have appreciable change with respect to the choice of t. The flavor-singlet potential is shown to 
have an "attractive core" and to be well localized in space. Because of the latter property, no significant 
volume dependence of the potential is observed within the statistical errors, as seen in Fig. [2TT Left). 
As the quark mass decreases in Fig. l^lTRight) , the long range part of the attraction tends to increase. 

The resultant potential is fitted by the following analytic function composed of an attractive Gaussian 

core plus a long range (Yukawa)^ attraction: V{r) = bie~^'^^ + 63(1 — e"''*'" ) f e~^^''/rj . With the five 
parameters, &i, 2,3,4,5, the lattice results can be fitted reasonably well with x^/dof ~ 1. The fitted result 
for L = 3.87 fm is shown by the dashed line in Fig. [2T](Left). 
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Figure 22: The bound state energy Eq = —B^ and the rms distance J {r"^) of the if-dibaryon obtained 
from the potential at L = 3.87 fm. (Left) {t — to) /a dependence at m^^ = 673 MeV. (Right) Quark 
mass dependence at (t — tQ)/a = 10. Taken from Ref. 



Solving the Schrodinger equation with the fitted potential in infinite volume, the energies and the 
wave functions are obtained at the present quark masses in the flavor SU(3) limit. It turns out that, 
at each quark mass, there is only one bound state with binding energy 30-40 MeV. In Fig. [22r Left). 
the energy and the root-mean-squared (rms) distance of the bound state are plotted in the case of 
{t — to)/a = 9, 10, 11 at rrips = 673 MeV and L = 3.87 fm, where errors are estimated by the jackknife 
method. Although the statistical error increases as t increases, we observe small changes of central 
values, which are considered as the systematic errors. Fig. l22r Right) shows the energy and the rms 
distance of the bound state at each quark mass obtained from the potential at L = 3.87 fm and 
(t — to) /a = 10. Despite the fact that the potential has quark mass dependence, the resultant binding 
energies of the if-dibaryon are insensitive in the present range of the quark masses. This is due to the 
fact that the increase of the attraction toward the lighter quark mass is partially compensated by the 
increase of the kinetic energy for the lighter baryon mass. It is noted that there appears no bound state 
for the potential of the 27-plet channel in the present range of the quark masses. 

The final results of the binding energy Bh and the rms distance 



(r^) are summarized below, where 
the 1st and 2nd parentheses correspond to statistical errors and systematic errors from the t-dependence, 
respectively. 

= 32.9(4.5)(6.6) MeV 

= 0.823(33) (40) fm 
= 37.4(4.4) (7.3) MeV 

= 0.855(29)(61) fm 
= 35.6(7.4)(4.0) MeV 

V2) = 1.011(63)(68) fm 

Since the binding energy is insensitive to the quark masses, there may be a possibility of weakly bound 
or resonant iJ-dibaryon even in the real world with lighter quark masses and the flavor SU(3) breaking. 
To make a definite conclusion on this point, the A A — NE — SS coupled channel analysis is necessary 
for H in the (2+l)-fiavor lattice QCD simulations, as discussed in Sec. 17.21 

Recently the existence of if-dibaryon is also suggested by a direct calculation of its binding energy in 
2-1-1 full QCD simulations [60], where Bh = 16.6(2.1)(4.6) MeV is reported in the L — ;■ oo extrapolation 
at m^ ~ 389 MeV . 
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8 Other applications 

In the last section, some other apphcations of the potential method are reviewed. 

8.1 Three nucleon force 

Recent precise calculations of few-nucleon systems clearly indicate that the 2 nucleon force alone is 
insufficient to understand the nuclei, which calls for three (and/or more) nucleon forces. Actually, the 
three nucleon force (TNF) is supposed to play an important and nontrivial role in various phenomena 
in nuclear and astrophysics. For the binding energies of light nuclei, the attractive TNF is required to 
reproduce the experimental data. On the other hand, the repulsive TNF is necessary to reproduce the 
empirical saturation density of symmetric nuclear matter. For the equation of state(EoS) of asymmetric 
nuclear matter, a repulsive TNF is required to explain the observed maximum neutron star mass. 

Pioneered by Fujita-Miyazawa [61], the TNF has been mainly studied from the two-pion exchange 
picture with the A-excitation. In addition, the repulsive TNF is often introduced phenomenologi- 
cally [62j. Recently, the TNF based on chiral effective field theory is developing [63J, but the unknown 
low-energy constants can be obtained only by the fitting to the experimental data. Since the TNF 
originates from the fact that the nucleon is not a fundamental particle, it is essential to study the TNF 
from fundamental degrees of freedom(DoF) , i.e., quarks and gluons. 

In Ref. [6l| 165] such first-principle calculations of the TNF in lattice QCD have been reported. 
If the potential method is applied to the three nucleon (3N) system, a straightforward calculation is 
impossible due to the significantly enlarged DoF. In Ref. [BH ES]) ^^o different approaches have been 
considered. 

The NBS wave function for 3N is defined by 

<^'^(n2,n23)e-'^* = {0\N{x,,t)N{x2,t)N{x;,t)\W) (120) 

where ri2 = Xi — X2, ri23 = x^ — {xi + X'^jl are the Jacobi coordinates, and \yV) is the 3N state with 
energy W . At the leading order of the velocity expansion, the NBS wave function satisfies 



"^m - TT V,\23 + Y. ^2N,ii(a;i - Xj) + yTNF(r'i2, r-ias) 



2/ii2 "' 2fi 



123 



«<J 



.W/-^ „ \ TP,„W 



V {ri2,ri23) = Eip (ri2,ri23) 



(121) 
where V2N,ii(a^i — Xj) denotes the potential between (z,j)-pair, Vtnf(^i2, ''"123) the TNF, fii2 = m7v/2, 
yUi23 = 2miv/3 the reduced masses. If {p^{ri2, r'123) is calculated for all r^, r^s and all V2N,ij(a5i — Xj) 
are available by lattice calculations, Vtnf(^i2, ''"123) can be extracted. Unfortunately, this is not the case: 
Since both ri2 and ri23 have L^ DoF, the calculation cost is more expensive by a factor of L^ compared 
to the 2N system. Furthermore, the number of diagrams to be calculated in the Wick contraction tends 
to diverge with a factor of A^„! x A^^! (A^^^^^ are numbers of u,d quarks in the system). It is also noted 
that not all 2N potentials are available in lattice QCD at this moment: Only parity-even 2N potentials 
have been obtained so far. 

The first method in Ref. [61] to avoid these problems is to consider the effective 2N potential in the 
3N system by taking the summation over the location of the spectator nucleon A^(a;3), 

^^(ri2) = E ^"^(^12, ^23) = E ¥^"^(^12, ^123). (122) 

3^3 '''123 

The effective potential between N{xi) and N{x2) is then defined by 



2/ii2 ■ ''" 



VL + Kff(n2) 



ip^{r^2) = Eip^{r^2)■ (123) 
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(Left) Effective 2N potentials, where red, blue, brown points correspond to V^gg 
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potential, respectively. (Right) The difference between the effective 2N and the 



genuine 2N for V^ 



/=o,s=i 



potential. Taken from Ref. [6 



In this calculation, the DoF of ri23 is integrated out beforehand, and thus the calculation cost is reduced 
by a factor of ~ l/L^, compared to the straightforward calculation. The difference Kff(^) — V2Ar(r) can 
be considered to be the "finite density effect" in the 3N system. Part of this effect is attributed to 
the genuine 2N potential with the nontrivial 3N correlation, while another originates from the genuine 
TNF. 

The triton channel( / = 1/2, J^ = 1/2+) is studied as the 3N system. Since the spectator nucleon 
is projected to the S-wave, the possible quantum numbers between the (effective) 2N are only ^'^"'"^Lj = 
^Sn, ^Si, ^Di. Gauge configurations in 2-flavor QCD on a 16^ x 32 lattice at a ~ 0.16 fm[66] are 



'0, 

employed for the calculation at m^ ~ 1.13 GeV and mi^ ~ 2.15 GeV. 

Fig. [23r Left) show results for Kfj(r) in the triton channel at t — to = 85 where the constant shift 
by energy is not included for the central potentials. It is noteworthy that \4flf(r) is obtained with good 
precision even though the signal to noise ratio is expected to be worse for more quarks in the system. 
Fig. [23r Right) gives V^flF(r) — V2n('") for the tensor potential, which is free from the constant shift. The 
difference is consistent with zero within a few MeV statistical errors. Similar results are reported for 
central potentials as well in Ref.[6l]. There is no indication of a TNF effect. A possible explanation is 
that the TNF effect is suppressed at heavy quark mass. Basically similar results are obtained, however, 
for lighter pion masses(m7r ~ 0.7 GeV and 0.57 GeV)[M]. Another possibility is that the TNF effect 
is suppressed by the summation over the location of the spectator nucleon. While the TNF effect is 
expected to be enhanced when all three nucleons are close to each other, such 3-dimensional spacial 
configurations have small contributions in the spectator summations. 

In order to assess this possibility, a second method has been investigated in Ref. [641 [65] . where the 
linear setup with ri23 = is used for the 3N wave function. In this case, the third nucleon is attached to 
(1, 2)-nucleon pair with only S-wave. Considering the total 3N quantum numbers of / = 1/2, J^ = 1/2+, 
the wave function can be completely spanned by only three bases, which can be labeled by the quantum 
numbers of (l,2)-pair as ^5*0, ^5*1, ^Di. Therefore, the Schrodinger equation is simplified to the 3x3 
coupled channel equations with the bases of v^i^g, V'ssi, ^P^Di- Even in this case the subtraction of V2N 
remains nontrivial: the parity-odd potentials, which must be subtracted, are not available in lattice 
QCD at this moment. The subtraction problem of parity-odd potentials can be avoided for triton by 
using the symmetric wave function. 



^s 



Vq 



p^n^ni + p^nin^ - n^riip^ + n^rz^p^ + n^p^rii - riip^n^ 



(124) 
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Figure 24: (Left) The wave function with hnear setup in the triton channel. Red, blue, brown points 
correspond to ips, fu, ^'-^d^i respectively. (Right) The scalar /isoscalar TNF in the triton channel, 
plotted against the distance r = |t'i2/2| in the linear setup. Taken from Ref. [6i 



Combined with the Pauli principle, it is automatically guaranteed that any 2N-pairs couple with even 
parity only, since this wave function is anti-symmetric in spin/isospin spaces for any 2N-pairs. Therefore 
the TNF can be extracted unambiguously in this channel, without the information of parity-odd 2N 
potentials. 

The same gauge configurations used for the effective 2N potential study are employed in the 
numerical simulations. Fig. l24r Left) gives each wave function of ips = 772 (""^^So + V'^sJ) fu = 
4=(+'i/'i5(, -|- ip3s-^), ipsjy-^ as a function of r = |r'i2/2| in the triton channel at t — to = 8. Among these 
three (ps dominates the wave function, since ips contains the component for which all three nucleons 
are in S-wave. 

By subtracting the V2N from the total potentials in the 3N system, the TNF is detemined. Fig. |2l] 
(Right) shows results for the scalar/isoscalar TNF, where the r-independent shift of V2N is not included, 
and thus about (9(10) MeV systematic error is understood. There are various physical implications in 
Fig. 121] (Right). At the long distance region of r, the TNF is small as is expected. At the short 
distance region, the indication of a repulsive TNF is observed. Recalling that the repulsive short-range 
TNF is phenomenologically required to explain the saturation density of nuclear matter, etc., this is 
a very encouraging result. Of course, further study is necessary to confirm this result, e.g., the study 
of the ground state saturation, the evaluation of the constant shift by energies, the examination of the 
discretization error. 



8.2 Meson-baryon interactions 

The potential method can be naturally extended to the meson-baryon systems and the meson-meson 
systems. In this subsection, two applications of the potential method to the meson-baryon system are 
discussed. 

The first application is the study of the KN interaction in the I{J^) = 0(l/2~) and l(l/2~) 
channels by the potential method. These channels may be relevant for the possible exotic state G^, 
whose existence is still controversial. 

The KN potentials in isospin 7 = and 1 = 1 channels have been calculated in 2 + 1 full 
QCD simulations, employing 700 gauge configurations on a 16^ x 32 lattice at a = 0.121(1) fm and 
{m^,mK,mN) = (871(1), 912(2), 1796(7)) in unit of MeVjSZ]. 

Fig. I25] shows the NBS wave functions of the KN scatterings in the 7 = (left) and 7=1 (right) 



43 



0.4 



-i i 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.1 

r[fm] 



0.8 



0.6 



— 0.4 



0.2 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.1 

r[fm] 



Figure 25: The NBS wave function of the KN scattering in the / = (left) and the 1 = 1 (right) 
channels. Taken from Ref. |67] . 
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Figure 26: The LO potential for the KN state without the energy shift E in the / 
1 = 1 channels (right). 



(left) and the 



channels. The large r behavior of the NBS wave functions in both channels do not show a sign of a 
bound state, though more detailed analysis is needed with larger volumes for a definite conclusion. On 
the other hand, the small r behavior of the NBS wave functions suggests a repulsive interaction at short 
distance (r < 0.3 fm). The repulsion in the 7 = 1 channel seems to be stronger than that in the 7 = 
channel. 

The LO potential V{r) for the KN state without the constant energy shift E is shown in Fig. [2Blin 
the 7 = (left) and 7 = 1 (right). As expected from the NBS wave functions in Fig. [23, the repulsive 
interactions are observed at short distance in both channels, while the attractive well appears at the 
medium distance (0.4 < r < 0.8 fm) in the 7 = channel. These results indicate that there are no 
bound states in 7(J^) = 0(l/2~) and 1(1/2") states at m^ ~ 870 MeV. 

The second application is to study the charmonium-nucleon interaction using the potential method. 
Since charmonium does not share the same quark flavor with the nucleon, the charmonium-nucleon 
interaction is mainly induced by the genuine QCD effect of multi-gluon exchanges. Theoretical studies 
based on QCD suggest that the cc-N interaction is weakly attractive. It is argued that the cc-nucleus (A) 
bound system may be realized for the mass number A > 3 ii the attraction between the charmonium 
and the nucleon is sufficiently strong [681 ES]- Precise information on the cc-N potential VccN{r) is 
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Figure 27: (Left) The effective central potential in the s-wave r]c-N system at m.^ = 640 MeV. The 
solid line is the fit by the Yukawa potential while the dashed line is the one by the phenomenological 
potential adopted in Ref. [68]. (Right) The volume dependence of the r]c-N potential. Taken from 
Ref. [71] • 



therefore indispensable for exploring nuclear-bound charmonium states such as the rjc-^He or J/ip-^He 
bound state in few body calculations [70]. 

In Ref. [72], the charmonium- nucleon potentials are calculated in quenched QCD on 16^ x 48 and 
32^ X 48 lattices at a ~ 0.94 fm at three different values of the light quark mass corresponding to 
(m^r, tun) ~ (640, 1430), (720, 1520), (870, 1700) in unit of MeV and one fixed value of the charm quark 



mass corresponding to m„ ~ 2920 MeV and m 



ij/* 



3000 MeV. 



The effective central rjc-N potential, evaluated from the NBS wave function with measured values 
of E and /x, is shown in Fig. [271(Left). The rjc-N potential clearly exhibits an entirely attractive 
interaction between charmonium and the nucleon without any repulsion at all distances. Absence of 
the short range repulsion (the repulsive core) may be related to absence of the Pauli exclusion between 
the heavy quarkonium and the light hadron. The interaction is exponentially screened in the long 
distance region r > 1 fm. The fit of the potential with the Yukawa form —76""^ /r gives 7 ~ 0.1 and 
a ~ 0.6 GeV (solid line in Fig. 1271) . which are compared with the phenomenological values 7 = 0.6 and 
a = 0.6 GeV in Ref. [HE] (dashed line). The strength of the Yukawa potential 7 is six times smaller 
than the phenomenological value, while the Yukawa screening parameter a obtained from lattice QCD 
is comparable with the phenomenological one. The cc-N potential is found to be rather weak in lattice 
QCD. 

As shown in Fig. [27l(Right). there is no significant difference between potentials at two different 
spatial sizes {La ~ 3.0 and 1.5 fm). The size dependence of the rjc-N potential seems small since the 
rjc-N potential is quickly screened to zero and turns out to be short ranged. 

No large quark- mass dependence is observed in Fig. [28](Left). This may be understood by the 
argument that the cc-N interaction is mainly governed by multi-gluon exchanges, which do not depend 
explicitly on the quark mass. Taking a closer look at the inset of Fig. [28l(Left). however, one finds 
that the attractive interaction in the r/c-N system tends to get slightly weaker as the light quark mass 
decreases. 

The spin- independent part of the central J/ip-N potential is shown at m.^ = 640 MeV in Fig.[28](Right). 
together with the r]c-N potential for comparison. While no qualitative difference between the t]c-N and 
J/ip-N potentials is observed, the attractive interaction in the J/ip-N potential is a little stronger than 
the r]c-N potential. The attraction, however, is still not strong enough to form a bound state in the 
J/ijj-N channel. 
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Figure 28: (Left) The quark-mass dependence of the ?7c-iV potentiaL (Right) The spin-independent part 
of the central J/ip-N potential at m^ = 640 MeV, together with the t]c-N potential for comparison. 
Taken from Ref. [7T] . 

8.3 Two-color QCD and potentials 

In Ref. [72], potentials between baryons in two-color (SU(2)) QCD are investigated on the lattice. 

In two-color QCD, two quarks (diquark) form a "baryon", whose local interpolating operator is 
given by 

Ai,r(^, t) = e'^'qiix, t)Tq,{x, t), (125) 

where e"''' is the 2x2 anti-symmetric tensor, i,j are flavor indices and T = C, C75, Cji^, 6*7^75 with the 
charge conjugation matrix C. The NBS wave function is then defined by 



¥'S(-)e-^* 



{0\T{D,,,rix + r,t)DkiAx,t)}\W) 



(126) 



where \W) is an eigenstate of two color QCD with total energy W and "baryon" number 2. The 
potential derived from this wave function is denoted as V^i^i{r). 

The potentials have been calculated in SU(2) quenched QCD on a 24^ x 48 lattice at a :ii 0.1 fm 
determined from the string tension with the assumption that ^/a = 440 MeV. The lightest "baryon" 
corresponds to the scalar diquark state (F = C). Therefore the NBS wave function is evaluated for 
F = C at four values of the quark mass, which give ms = 1.044(2), 0.836(2), 0.618(2), 0.377(3) in lattice 
unit with ms being the scalar diquark mass. The potential is then extracted from the Ai state with 
T = C, and is therefore denoted by Vij^ki{r). 

Fig. [2^ shows the LO central potentials plotted as functions of r for {i,j,k,l) = (1, 2, 3, 4) (left) 
and (1, 2, 1, 2) (right), where quark-exchange diagrams are absent for the former while they exist for the 
latter . 

The potential Vi2,34('") has an attractive interaction at all length scales, which becomes stronger as 
ms decreases. The r-dependence is monotonic at all values of m,s- At large r (r > 4 in lattice units), 
the potential has small m,s dependence, while it has stronger m,s dependence at small r (r < 4). The 
magnitude of the potential decreases as m^s increases and finally the potentials at ms = 0.836 and 1.044 
coincide with each other. 

For Vi2,i2, on the other hand, strong repulsions appear at short distance. The m^s dependence is not 
monotonic: The potential is a smooth function of r at small m^, while it has a pocket at intermediate 
distance for m^ = 1.044. As m^s decreases, the repulsive core rapidly increases and the attractive pocket 
disappears. 

The qualitative difference between the two cases suggests that contributions from quark exchange 
diagrams, which more or less represent the Pauli exclusion effect, are responsible for the repulsive core. 
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Figure 29: The LO central potentials Vi2,34{r) (Left) and Vi2,i2(^) (Right) as functions of the relative 
distance r in lattice unit. Here k = 0.135, 14, 145 and 0.15 correspond to ms = 1.044,0.836,0.618 and 
0.377 in lattice unit, respectively. Taken from Ref. [72]. 



On the other hand, the attraction may be explained by gluon exchanges, which are the main part of 
the potential V^i2,34('")- More details analysis can be found in Ref. 



9 Conclusion 

In this report, we have reviewed the progress on the study of hadron interactions via the potential 
method in lattice QCD. The key quantity of the method is the NBS wave function of the two particle 
state. Not only is the asymptotic behavior of the NBS wave related to the scattering phase shift 
(the phase of the S-matrix) in QCD, but also the non-local but energy-independent potential (or the 
interaction kernel) can be extracted from the NBS wave function in the non-asymptotic region. By 
construction the potential correctly reproduces the scattering phase shift at all energies below the 
inelastic threshold. In practice the non-local potential is approximated by the velocity expansion in 
terms of local functions, so that physical observables such as the scattering phase shift are approximately 
calculated. It is also possible to check the accuracy of the approximation. 

Using the local potential approximation at the leading order in the velocity expansion, nucleon- 
nucleon, hyperon-nucleon and hyperon-hyperon interactions have already been investigated successfully. 
The BB potential in the flavor SU(3) limit shows the attraction in the flavor singlet channel, which is 
strong enough to form a bound state, if-dibaryon in this limit. The repulsive core is analyzed in terms 
of the operator product expansion and the renormalization group. 

The potential method is extended to the case that the total energy is above the inelastic threshold. 
This method seems to work for the BB interaction with S = —2 and / = 0, where the coupled channel 
potentials among AA, NS and SS are obtained. The potential method is also applied to the three 
nucleon force, the meson-baryon interactions and two color QCD. 

Of course, more careful studies of systematic errors such as the finite volume effect, dynamical quark 
effect, quark mass dependence and the lattice spacing effect are needed, before applying potentials 
obtained in lattice QCD in order to investigations of nuclei and hyper-nuclei. 

Finally I stress here that the potential derived from the NBS wave function adds a new tool to inves- 
tigate hadron interactions in lattice QCD, in addition to the standard finite volume method proposed 
by Liischer [14j. In particular, the hadron scattering above the inelastic threshold can be treated in 
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lattice QCD by the potential method. Other extensions of this method will also be looked for. 
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A Asymptotic behaviors of the NBS wave function in the spin 
triplet channel 

There are three components in the spin triplet channel at given J . One component corresponds to the 
/ = J scattering, which is extracted by taking \s1S2) = 75 (l + 2' ^"2) ~ I ~ 2' ~\n- As r — )■ 00, this 
component becomes 






(127) 



where 



Z'''{S) = ^DUn^)UiV)Ui-V)J2xiS,s,){l,S,h,s,\l,h-s,) 

X (^(/, h - s„ /, S\l, h - s„ +-, +-) - (/, /. - s„ I, S\l, h - s„ --, --)) . (128) 



The / = J=f1 components, which correspond to \s1S2) = 4j (||, — |) i I ~ |i +2)) = 1='=)' ^^^ niixed 
with each other, and those asymptotic behavior is given by 



J,M,l,lz 



5i,j-i 
Sij+i 



^Ri'US = hk) 



(129) 



where 



Z-^^iS) = ZUiy)U{-W)Y.x{S,s,){lShs,\JM) (130) 

Sz 

Lif{S,k)t, = E {JMlS\JMs,S2)[X0{k)],^_,^^,^ (131) 

Sl — S2=±l 

Ri^.s,iS,kh = [0-\k)x] Di^^^^JQ^), (132) 



and Si — S2 becomes ±1 in the above equation. Here ti = ±1 corresponds to |±) introduced before, and 

1 _i\ |_1 +1^ 

2' 2' ' I 2' ^2/ 



X, which transforms the helicity basis {| + |, —k), I — |) +|)} to the new basis {|+), |— ), is given by 



A- = A- = -L(; _lj. (133) 
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The mixing matrix in the new basis has the form 

0{k) = 
where Oj{k) is the mixing angle between / = J ± 1 



cos 9 j{k) —sin 9 j{k) 
sm9j{k) cos 9 j{k) 



(134) 



B Octet baryons 

Local octet baryon operators are defined by 

na{x) = eabc{u"-{x)C-frjd^{x))dl{x) 

-- -eabc{u''{x)C-f5s\x))Ua{x) 



S+(:r) 



-eabc-^[{d''{x)C^,s\x))ul{x) + u''{x)C^5s\x))d:,{x)] 



S-(x) = -eabcid^ix)C^,s\x))dl{x) 
S°(x) = eabcis''ix)C^,u\x))slix) 
E-{x) = eabc{s''{x)C-f5d''{x))s''^{x) 

Aa{x) = ~eabc^[{d''{x)C^,s\x)X^{x)+u-{x)C^,s'{x))d'^{x)-2u^{x)C^,d\x))s'^^^^^ 

from which the irreducible representations of the flavor SU(3) with baryon number B = 2 are constructed 
as 



55) (27) 


bbY^^^ 


{BBf^^ 


BB)'^^ 


55)(io) 


BBY^'^^ 



-AA- 



'iAA 



-AA 



-SS 



'5ee 



-SS 



'-NE 



•-pn 



2^" 



-np 



1^- 



-p 






1 „ 



-n, 



■^0 



'U, 



■^0 



n 



(136) 



where SS and NE are defined 



NE 



-S+S" 



4^" 



Ov^O 



-S^S' 



i^ P 



-S"S" 



4 



'iE«„. 



(137) 



Unitary matrices which rotate the fiavor basis to the baryon basis are given as follows. 
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1. S=— 1, 1=1/2, spin-singlet. 




2. S=-l, 1=1/2, spin-triplet. 




_9_ 
10 
J_ 
10 



3. S=— 2, 1=0, spin-singlet. 




V 



/27 
40 
1 



40 
40 



4. S=— 2, 1=1, spin-singlet. 




5. S=— 2, 1=1, spin-triplet. 




V 



6. S=— 3, 1=1/2, spin-singlet. 




7. S=— 3, 1=1/2, spin-triplet. 
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